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ABSTRACT 


Many  applications  of  wind-wave  forecasting  require 
a completely  objective,  high-speed  method  of  fore- 
casting wave  spectra  for  large  areas  simultaneously. 
Consequently,  a numerical  model  for  forecasting  wave 
spectral  fields  was  defined,  programmed,  and  tested 
in  two  sample  situations. 


The  model  uses  the  Neumann  spectrum  with  angular 
dispersion  as  defined  by  the  Stereo  Wave  Observa- 
tion Project.  The  usual,  problems  in  defining  fetch 
shapes  and  motions  are  eliminated  by  consideration 
of  an  equivalent  Lagrangian  duration. 


In  both  test  computations  the  method  adequately 
forecasts  the  rise  in  significant  height  as  the 
winds  increase.  However,  the  spectral  frequencies 
forecast  are  too  high  and  the  swell  does  not  decay 
fast  enough,  k. 
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FOREWORD 


This  research  hegan  under  the  sponsorship  of  the  Office  of  Naval  Research, 
Contract  NOnr  285(03)  In  the  Department  of  Meteorology  and  Oceanography, 

New  York  University.  Some  machine  time  and  programming  assistance  was  pro- 
vided hy  the  Numerical  Weather  Prediction  Unit  of  U.  S.  Navy  Fleet  Weather 
Central  which  has  since  become  the  Fleet  Numerical  Weather  Facility  at 
Monterey,  California.  Machine  time  for  the  final  programming  and  running 
was  provided  by  the  Lockheed  Missiles  and  Space  Company,  a Division  of 
Lockheed  Aircraft  Corporation,  under  Contract  Nord  17017*  'All  of  these 
organizations,  and  people  responsible  for  machine  operation  and -programming, 
deserve  recognition  and  appreciation. 

Special  recognition  and  appreciation  are  due  Professor  W.  J.  Pierson,  Jr., 
who  has  encouraged  this  effort  from  its  very  inception  and  who  has  helped 
immeasurably  through  many  discussions,  letters,  and  criticisms.  Appreciation 
is  also  accorded  to  Messers.  Lionel  Moskowitz  and  Masashi  Murakami  who  scaled 
the  winds  from  the  synoptic  charts  for  Case  II. 

This  report  is  essentially  a copy  of  a thesis  submitted  to  the  Department  of 
Meteorology  and  Oceanography  New  York  University. 
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Chapter  1 
INTRODUCTION 


1-1  IMPORTANCE  OF  AND  NEED  FOR  OCEAN  WAVE  FORECASTING 


Accurate  large-scale  ocean  wave  forecasting  is  extremely  important.  High  seas 
can  damage  large  ships  and  be  a real  danger  to  smaller  ones.  Seaplanes  can 
only  land  in  relatively  low  seas.  Hydrofoils  and  surface  effect  machines 
(ground  effect  machine)  cannot  operate  in  rough  seas.  Waves  often  force  ships 
to  change  speed  and  heading.  The  design  of  ships  and  shore  facilities,  such 
as  sea  walls  and  docks,  must  consider  the  forces  that  are  exerted  by  the  ex- 
treme waves.  Wave  forecasting  can  allow  ship  operators  to  plan  ahead,  thus 
saving  time  and  money,  while  also  decreasing  danger,  damage,  and  discomfort. 


There  are  many  military  applications  in  which  wave  forecasts  are  useful.  Be- 
sides needing  forecasts  for  problems  that  are  similar  to  the  commercial  uses, 
there  are  many  other  needs.  For  example,  rendezvous  for  refueling  or  trans- 
ferring of  men  or  equipment  have  a double  need  for  wave  forecasts.  The  speed 
of  the  ship  toward  the  rendezvous  is  affected  by  sea  conditions,  then,  at  the 
rendezvous,  high  seas  can  hinder  or  halt  the  transfer.  Waves  are  of  critical 
importance  to  aircraft  carrier  operations  because  high  seas  cause  large  ship 
motions  in  which  airplanes  cannot  land  safely.  Surface  waves  are  important 
to  submarines  and  anti-submarine  crews  because  the  noise  generated  by  the  waves 
can  provide  a perfect  hiding  place.  Such  things  as  antennas  and  snorkels  can 
be  hidden  in  the  clutter  on  radar  screens.  The  sea  condition  is  one  of  the 
most  important  and  critical  environmental  factors  in  almost  any  naval  or  ma- 
rine operation. 


Since  it  is  difficult  to  make  measurements  of  ocean  waves,  hindcasting  tech- 
niques allow  computation  of  wave  statistics  for  design  purposes.  In  order  to 
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solve  these  and  other  practical  problems,  the  entire  spectrum  of  the  sea  and 
swell  must  be  known  and  used.  The  reasons  for  this  are  now  fairly  well  ac- 
cepted and  are  also  explained  from  a practical  viewpoint  by  D.M.  Aspinwall 
(Reference  1,  i960),  among  others. 

Another  more  important,  more  basic,  and  longer-range  problem  exists.  This  is 
the  problem  of  advancing  our  knowledge  of  how  the  wind  transfers  energy  to  the 
water  and  how  different  wave  trains  interact  to  transfer  and  dissipate  energy. 
This  problem  has  been  studied  most  recently  by  O.M.  Phillips  (Reference  2,  1961), 
L.J.  Tick  (Reference  3,  1961),  W.J.  Pierson  (Reference  k,  1961),  O.M.  Phillips 
and  E.J.  Katz  (Reference  5,  1961),  and  K.  Hasselmann  (Reference  6,  i960),  among 
others,  frcm  the  theoretical  viewpoint.  Many  others  have  attacked  the  air-sea 
boundary  layer  problems,  such  as  wind  stress.  However,  no  one  has  yet  estab- 
lished an  empirical  model  by  which  this  energy  exchange  can  be  studied  objec- 
tively. Only  by  the  establishment  of  objective  models,  such  as  the  one  pro- 
posed here,  can  the  future  theoretical  approaches  be  verified  (as  yet  none  has 
reached  the  stage  of  readiness  for  verification).  It  is  also  hoped  that  this 
model  will  help  with  such  theoretical  problems  as: 

1.  The  interaction  of  opposing  wave  trains. 

2.  The  proper  fetch  or  duration  needed  for  "full  development. " 

3.  The  true  shape  (or  shapes)  of  the  spectrum. 

4.  The  effect  of  wind  on  swell. 

5.  The  importance  of  non-linear  effects  such  as  breakers. 

Some  of  the  more  important  practical  problems  in  which  a model,  such  as  the 
one  proposed,  would  be  helpful: 

1.  Hindcasting  of  wave  spectral  statistics. 

2.  Hindcasting  to  learn  the  area  distribution  characteristics,  i.e.,  the 
joint  probability  of  sea  conditions  occurring  with  different  intensities 
at  several  locations . 
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3.  Forecasting  for  ship  routing. 

4.  Forecasting  for  naval  development. 

5.  Improving  present  graphical  methods  of  wave  forecasting. 

Verification  of  a model  such  as  this  cannot  he  complete  until  two-dimensional 
spectra  are  available  for  waves  measured  on  a synoptic  basis.  In  the  mean- 
time, rough  verification  has  been  carried  out,  first,  with  the  estimated  sig- 
nificant heights  from  synoptic  observations  and,  second,  with  the  spectra  of 
a particularly  severe  sea  condition  recorded  by  the  O.W. S.  Weather  Reporter. 
If  the  National  Institute  of  Oceanography  in  Great  Britain  had  not  developed 
ship-borne  wave  recorders  (M.J.  Tucker,  Reference  7,  1956)  and  used  them  for 
systematic  synoptic  observations,  even  this  much  verification  could  not  have 
been  attempted. 


1-2  PRESENT  TECHNIQUES 

Most  of  the  important  graphical  methods  of  forecasting  wave  spectra  were  sum- 
marized and  discussed  by  G.  Neumann  and  W.J.  Pierson,  Jr.  (Reference  8,  1957)* 
Since  that  time,  J.  Derbyshire  (Reference  9,  1959)  has  modified  his  method, 
and  the  C.L.  Bretschneider  (Reference  10,  1959)  and  the  "Densites  Spectro- 
Angulaires"  (DSA)  of  R.  Gelci,  H.  Cazale,  J.  Vassal  (Reference  11,  1957)  methods 
have  been  proposed.  W.J.  Pierson,  Jr.  (Reference  12,  1959),  H.  Walden  (Ref- 
erence 13,  1961),  and  others  have  compared  several  forecasting  methods  with 
inconclusive  results.  (See  also  W.J.  Pierson,  Jr.  and  G.  Neumann,  Reference  14, 
1961).  With  the  exception  of  a recent  application  by  J.  Darbyshire  (Ref- 
erence 15,  1961)  and  a third  model  of  DSA  ( R.  Gelci  and  P.  Chavy,  Reference  16, 
1961),  all  of  these  methods  are  useful  only  for  a limited  number  of  points. 
Darbyshire,  in  opposition  to  others,  maintains  that  duration  and  fetch  are 
relatively  unimportant.  He  can,  therefore,  easily  use  his  method  for  comput- 
ing a two-dimensional  map  of  wave  conditions.  This  work  by  Darbyshire  was 
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carried  out  while  the  present  study  was  in  its  final  stages.  The  DSA  numerical 
method,  which  is  presently  being  studied  in  France,  is  quite  similar  in  princi- 
ple to  the  methods  described  herein,  except  for  the  simplification  that  the  en- 
ergy increment  for  each  frequency  component  which  is  added  does  not  depend  on 
the  complete  initial  spectrum,  but  only  on  the  initial  energy  within  the  same 
component. 

\ J 

There  has  been  a great  deal  of  controversy  as  to  which  of  the  proposed  theoret- 
ical spectra  is  best.  No  one  method  has  been  definitely  proven  to  be  best. 
Since  the  W.J.  Pierson,  Jr.  - G.  Neumann  - R.W.  James  (Reference  17,  1955) 
method  is  widely  accepted  in  this  country,  it  was  used  in  the  present  study. 


The  Pierson-Neumann-James  system  and  most  other  wave  forecasting  methods  use 
three  basic  parameters  to  forecast  the  sea  conditions.  These  are;  the  wind 
speed,  the  area  over  which  the  high  winds  are  fairly  constant,  called  "fetch", 
and  the  time  duration  that  the  wind  blows.  Earlier  methods  did  not  consider 
the  width  of  the  fetch  important;  however,  Pierson,  Neumann,  and  James  showed 
a necessity  for  measuring  this  width  to  define  the  swell  filters.  These  fil- 
ters are  a simple  way  of  keeping  up  with  the  energy  content  in  the  various 
frequency-directional  increments  in  the  spectrum.  The  Pierson-Neumann-James 
system  was  extended  by  L. S.  Simpson  (Reference  18,  1955)  to  allow  for  some 
movements  of  the  fetch.  He  showed  large  effects  on  wave  height  caused  by  the 
movement  of  the  fetches;  however,  the  method  is  laborious  and  considers  only 
simple  motions  of  the  fetch. 

Since  there  is  a great  demand  for  a two-dimensional  map  of  ocean  waves,  the 
techniques  of  Pierson-Neumann-James  have  been  used  by  the  United  States  Navy 
Hydrographic  office  and  by  others  to  forecast  re a waves  for  many  points  from 
which  plotted  map  contours  can  be  drawn.  However,  in  doing  this  work  by  hand, 
the  volume  of  the  work  usually  forces  the  analyst  to  neglect  swell  and  other 
fine  points.  W.E.  Hubert  (Reference  19,  1957)  has  used  an  IBM  70^;  the  co- 
cumulative  spectra  of  Pierson,  Neumann  and  James,  and  a large  number  of  grid 
points,  to  make  a forecast  map  of  significant  heights.  He  had  also  been  forced 
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to  neglect  swell.  Neither  of  these  mapping  methods  considers  moving  fetches 
of  changing  wind  fields. 

Another  difficult  problem  Is  In  delineating  a fetch.  In  most  actual  cases  a 
fetch  does  not  have  clear-cut  boundaries.  The  wind  Is  not  really  constant 
within  the  fetch.  Also,  there  Is  some  finite  transition  zone  around  the  edges. 
Two  Independent  analysts  will  seldom  define  the  same  fetch.  Therefore,  In  de- 
vising a method  to  make  two-dimensional  wave  forecasts,  It  would  be  convenient 
if  the  measurement  of  fetch  could  be  eliminated. 
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Chapter  2 
PHYSICAL  MODEL 


2-1  ASSUMPTIONS 

The  method  of  forecasting  used  Is  based  on  certain  physical  assumptions  that 

are  similar  to  those  used  by  Pierson,  Neumann,  and  James.  These  assumptions 

are  listed,  with  the  differences  from  Pierson,  Neumann,  and  James  underlined, 

as  follows: 

1.  The  wind  generates  the  wave  spectrum  from  high  to  low  frequency.  If  there 
Is  a gap  In  the  spectrum,  this  gap  must  be  completely  filled,  starting  at 
the  high  frequency  end  of  the  gap,  before  the  rest  of  the  spectrum  will  be 
developed. 

2.  The  wind  generates  component  waves  according  to  Neumann's  wave-generation 
graphs  with  angular  dispersion  as  given  by  the  results  of  Project  SWOP, 
L.J.  Cote  et  al,  (Reference  20,  i960). 

3.  The  annemometer-height  wind  velocity  at  a point  In  time  and  space  Is  re- 
presentative of  the  wind  within  a finite  area  surrounding  the  point  In 
space  and  time.  Further,  In  a practical  application,  this  wind  can  be 
found  with  sufficient  accuracy. 

4.  Wave  conditions  are  approximately  constant  over  a finite  area  around  each 
grldpoint. 

5.  The  linear  theory  of  water  waves  in  Infinitely  deep  water  holds,  so  that 
the  speed  of  propagation  of  the  component  waves  is  a function  of  only  the 
frequencies  of  the  components.  This  implies  a "linear"  spectrum. 

6.  There  Is  no  dissipation  of  wave  energy  In  deep  water,  also,  there  Is  no 
transfer  of  energy  between  components.  Thus,  for  example,  cross  seas  do 
not  Interact.  These  assumptions  are  not  necessary  to  the  development,  but 
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follow  Neumann  because  no  relationships  have  yet  been  established  that 
can  be  used  with  confidence. 


7. 


The  total  energy  increase  In  a unit  time  Interval  is  a function  of  (a)  the 
initial  energy  of  the  waves  already  moving  within  +90  degrees  to  the  wind 


direction  and  (b)  the  wind  speed.  This  particular  assumption  is  also  not 
necessary  to  the  development  and  could  be  revised  when  the  relationship  is 
better  understood. 


8. 


Wave  components,  once  generated,  move  along  approximate  great  circle  paths 
without  spreading  or  changing  direction. 


Although  turbulent  pressure  and  momentum  fluctuations  in  the  air  next  to 
the  surface  are  thought  to  be  important,  these  and  any  other  Important 
parameters  are  assumed  to  be  directly  proportional  to  the  "surface  wind 
velocity"  always. 


2-2  DESCRIPTION 


The  system  developed  in  this  paper  permits  one  to  dispense  with  an  Eulerian 
Fetch  if  the  duration  is  measured  in  "Lagrangian”  coordinates.  The  tern 
"Lagrangian"  is  used  here  to  mean  that  the  coordinates  used  in  the  forecast 
move  with  the  group  velocity  of  the  component  waves.  The  "Lagrangian"  dura- 
tion is,  therefore,  the  length  of  time  that  the  wind  is  adding  energy  to  a 
particular  moving  energy  packet  or  wave  group.  In  other  words,  the  distance 
that  the  coordinates  themselves  travel  while  the  waves  are  growing  is  equal 
to  the  fetch  distance  if  the  fetch  does  not  move.  This  is,  therefore,  a gen- 
eralization of  Simpson's  (Reference  18)  analysis  for  a moving  fetch. 


Thus,  the  speed  with  which  the  waves  will  move  through  the  fetch  actually 
limits  the  time  that  the  wind  can  add  energy  to  the  waves.  The  same  effect 
can  be  shown  in  a rough  way  by  using  the  co-cumulative  spectra  (CCS),  curves, 
developed  by  Neumann,  and  the  integrated  group  velocity  from  classical  linear 
wave  theory:  R = t/0.66f.  Where  R is  the  distance  of  travel  in  nautical  miles, 
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t is  time  in  hours,  and  f is  frequency  (reciprocal  period)  in  sec.-^.  To  show 
this,  one  can  calculate  the  distance  that  the  leading  edge  of  the  energy  will 
move  in  a finite  time  from  Neumann's  CCS  curves  as  shown  in  Table  2-1.  In  this 
table,  a very  rough  time-weighted  average  of  the  minimum  frequencies  was  pre- 
pared by  averaging  the  values  of  f at  five  points  equally  spaced  in  time. 

This  has  been  used  to  compute  an  equivalent  fetch  which  is  quite  close  to  the 
empirical  value  for  the  fetch.  Neumann  went  through  an  equivalent  process  in 
establishing  his  curves,  and  O.M.  Phillips  (Reference  21,  1958),  has  described 
the  process  theoretically. 


Table  2-1  Comparison  of  Lagrangian  Travel  Distance 
With  Eulerian  Fetch 


Wind  Speed  (kts) 

30 

40 

50 

Minimum  Duration  (hrs) 

23 

42 

69 

*/  -1\ 

Time-meaned  Frequency  (sec  ) 

0.1270 

0.0857 

0.0705 

Distance  Traveled  (nm) 

274 

742 

1483 

Minimum  Fetch  From  CCS  Curves 
(nm) 

280 

711 

1420 

jh  Error 

2.1 

4.4 

4.4 

♦Average  of  five  points  equally  spaced  in  time. 

It  should  be  stressed  that  no  assumption  was  made  regarding  the  shape  or  con- 
sistency of  any  "fetch".  This  is  because  in  this  model  the  actual  irregular 
and  changing  wind  field  as  a function  of  time  is  used.  The  method  of  specifi- 
cation alone  limits  the  area  and  temporal  irregularities  allowed. 

By  considering  the  above,  a method  for  allowing  the  wave  groups  to  move  through 
the  wind  field  at  their  own  propagation  velocities  has  been  devised.  Thus,  the 
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Chapter  3 
COMPUTING  SCHEME 

3-1  COMPUTING  PHILOSOPHY  AND  EQUIPMENT 

The  section  describing  the  physical  model  was  purposely  kept  very  brief  because 
the  usefulness  and  complexity  of  models  of  this  nature  must  depend  on  available 
computing  facilities.  This  problem  was  first  attempted  in  1958  with  only  a set 
of  tables  describing  the  spectrum  and  the  intention  of  using  graphical  addition 
on  a series  of  charts.  Each  chart  was  to  represent  a particular  spectral  com- 
ponent. This  might  have  answered  a few  questions  about  the  model,  but  could 
not  have  helped  prepare  for  future  advances.  The  volume  and  length  of  the  work 
was  already  near  the  maximum  that  a person  can  be  expected  to  perform.  Thus, 
it  was  decided  to  program  the  problem  for  a large-scale  electronic  computer. 

The  model  has,  therefore,  been  modified  slightly  by  characteristics  of  the 
available  computer.  The  scheme  described  here  is  that  used  on  an  IBM  7090 
having  32,768  words  of  core  storage.  The  problem  was  also  programmed  for  an 
IBM  704,  having  8,192  words  of  core  storage  plus  drums,  and  a CDC  1604.  Each 
was  slightly  different. 


3-2  COORDINATES 

The  coordinate  system  must  be  described  before  proceeding  further.  A grid  for 
the  area  spacing  and  individual  time  steps  for  the  temporal  coordinate  was 
chosen  as  approximately  square  with  2 degrees  of  latitude  on  a side,  and  a 
time  increment  of  two  hours.  The  grid  used  is  shown  in  Figure  3-1*  These 
choices  are  consistent,  since  the  motion  of  the  fastest  wave  component  is  such 
that  it  will  not  travel  to  the  next  gridpoint.  Implicit  in  the  choice  of  this 
coordinate  system  is  that  forecasts  will  not  be  more  accurate  in  space  than 
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the  gridpoint  spacing.  Ways  to  Improve  this  were  considered  too  complicated 
for  a first  model,  hut  are  discussed  at  the  end  of  this  paper.  This  spacing 
was  also  chosen  to  be  consistent  with  the  accuracy  of  wind  reports.  Wind  re- 
ports, in  general,  are  available  only  at  six  or  twelve-hour  increments;  and 
the  location  of  the  reporting  ship,  allowing  for  an  equivalent  distance  error 
Instead  of  observational  error,  Is  normally  less  than  two  degrees  latitude. 

The  location  of  Important  phenomena,  6uch  as  cold  fronts,  however,  is  often 
no  better  than  two  degrees  latitude.  Therefore,  a compromise  on  two  degrees 
and  two  hours  is  reasonable. 

It  was  stated  that  the  chosen  grid  was  square  and  approximately  two  degrees 
of  latitude  on  a side.  The  reasons  for  approximating  deserve  further  dis- 
cussion. Since  waves  travel  along  a great  circle  path,  a Gncmonic  projection 
would  probably  be  useful.  Since  this  projection  is  not  equal  in  area,  loca- 
tion effects  must  be  considered,  causing  a further  complexity.  Wave  trains 
leaving  from  one  gridpoint  at  a direction  of,  say,  30  degrees  will  arrive  at 
another  gridpoint  at  perhaps  43  degrees.  Thus,  the  computing  system  must  keep 
up  with  all  these  variations.  Because  the  model,  as  programmed,  uses  almost 
all  of  the  capacity  of  the  available  computers,  because  this  is  only  a first 
approximation  and  because  the  present  program  considers  only  the  middle  lati- 
tudes, a compromise  on  30  - 60  Lambert  Conformal  mapping  was  made.  On  this 
projection,  great  circles  are  not  quite  straight,  but  for  the  purposes  of  this 
study,  they  can  be  treated  as  straight.  In  middle  latitudes  a location  accuracy 
pf  only  +2  degrees  latitude  is  claimed.  The  Lambert  projection  also  has  a 
variable  spacing  of  the  latitudinal  grid,  but  this  is  only  a minor  effect. 

The  grid  used  is  treated  as  if  all  squares  were  exactly  the  same  size.  All 
directions  are  referred  to  the  coordinates  of  the  grid  rather  than  to  the 
coordinates  of  the  Earth.  Thus,  the  waves  follow  approximately  a great  circle 
by  moving  in  a constant  direction  on  the  grid.  No  consideration  was  given  to 
the  divergence  caused  by  the  curvature  of  the  Earth.  Ways  to  improve  these 
assumptions  are  considered  at  the  end  of  this  report. 
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3-3  SPECIFICATION  OF  THE  SPECTRUM 

The  second  decision  that  must  be  made  lies  in  specifying  the  spectrum.  There 
are  at  least  three  possible  ways  to  specify  the  spectrum  in  the  memory  of  the 
machine.  The  first  would  be  to  remember  all  of  the  generating  conditions  in 
such  a way  that  they  would  be  carried  along  with  the  wave  trains.  This  is, 
effectively,  the  method  used  for  graphical  computations.  However,  it  seemed 
cumbersome  when  considering  over  500  grid  points  simultaneously.  A second 
approach  might  be  to  fit  seme  polynomial  to  the  spectrum  at  each  gridpoint 
which  would  be  changed  as  the  wave  conditions  change.  This  approach  appears 
much  better  but  was  not  used  because  it  is  relatively  complicated  in  compari- 
son to  the  third  approach.  It  should,  certainly,  be  considered  for  future  at- 
tempts. Instead,  a simple,  straightforward  technique  was  used.  Ten  frequency- 
increments  and  twelve  directional-increments  were  chosen.  The  numerical  value 
of  the  energy  within  each  of  these  120  increments  is  remembered.  The  midpoints 
of  the  frequency  intervals  were  finally  chosen  as  having  equidistantly  spaced 
periods  of  2.85,  5.31,  7.77,  10.23,  12.69,  15.15,  17.61,  20.07,  22.53,  and 
25.00  seconds.  The  reasons  for  using  constant  period  increments  are  explained 
later.  Several  other  values  were  tried,  but  the  above  increments  gave  the 
best  results  without  using  refined  integration  techniques.  In  accordance  with 
Pierson,  Neumann,  and  James  (Reference  17)  and  O.M.  Phillips  (Reference  27, 
1961),  it  is  especially  important  to  have  very  high  frequencies  specified. 

This  is  because  energy  is  added  much  more  slowly  to  very  small  high  frequency 
waves  than  to  larger  waves.  Since  these  high  frequency  waves  propagate  so 
slowly,  and  because  there  is  no  dissipative  force  in  the  present  model,  the 
model  ocean  can  be  expected  to  become  filled  with  high  frequency  wave  com- 
ponents . 

The  directions  were  chosen  in  30  degree  Increments  (+15  degrees)  because  it  was 
believed  that  the  wave  directions  should  be  somewhat  less  accurate  than  the 
directions  of  the  winds  causing  them.  Surface  winds  are  normally  specified 
to  16  points  of  the  compass  (22.5  degrees)  or  to  10  degrees. 
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As  stated  above,  discrete  increments  have  been  chosen  to  represent  the  spectrum. 
This  causes  certain  errors  that  must  be  considered.  After  many  time  steps, 
swell  will  be  separated  so  that  no  energy  is  at  some  Intermediate  grid  point. 

In  other  words,  the  longer,  faster  waves  will  outrun  the  waves  of  the  next 
discrete  frequency  band  by  more  than  the  distance  to  the  next  grid  point.  For 
example,  in  the  extreme  case  with  a spectrum  having  energy  at  f = 0.04  and 
f = 0.0 6 the  increments  would  separate  at  the  rate  of  0.19  spaces  per  time 
step  and  would  be  two  spaces  apart  after  about  eleven  time  steps.  Similarly, 
the  discrete  directions  will  keep  energy  from  reaching  intermediate  grid  points 
after  moving  only  to  the  third  diagonal  grid  point,  which  could  be  as  rapid  as 
six  time  steps.  Even  though  these  errors  may  be  significant  in  same  cases, 
they  are  not  nearly  as  important  in  the  practical  case  because  fetches  are 
normally  larger  than  one  grid  point  in  length  and  width.  With  longer  fetches 
these  vacant  intermediate  grid  points  are  filled  from  adjoining  areas  of  the 
fetch.  Since  this  does  not  allow  the  energy  to  spread  out  but  keeps  it  con- 
centrated along  a given  path  and  at  a particular  frequency,  the  model  is  biased 
toward  higher  swells.  In  the  case  of  small  classical  fetches,  this  bias  would 
cause  the  forecast  field  to  be  irregular  with  adjacent  grid  points  being  either 
too  high  or  too  low. 

Specification  of  a map  of  the  final  forecast  spectra  is  harder  than  that  of 
the  initial  spectra.  To  present  120  values  at  each  of  519  grid  points  would 
confuse  the  reader  through  mere  volume.  Area  effects  would  not  be  apparent. 

It  Is,  however,  a simple  matter  to  print  any  parameters  of  interest  because 
all  values  are  in  the  machine.  It  would  be  easy  to  present  maps  of  total 
energy,  significant  height,  most  important  direction,  lowest  frequency,  fre- 
quency having  the  most  energy,  frequency  and  height  of  the  largest  swells, 
etc.  The  problem  is  presently,  there  is  little  data  with  which  to  verify  these 
maps,  and  no  use  for  most  of  them.  Therefore,  only  significant  height  is  pre- 
sented in  map  form  in  this  report.  The  complete  spectra  for  selected  grid 
points  are  printed  separately.  More  detail  on  this  is  given  in  the  section 
on  results  of  this  report. 
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3-4  PROPAGATION 

This  model  depends  upon  allowing  the  wave  energy  to  he  propagated.  Originally, 
a system  similar  to  the  Jacobian  method,  used  In  numerical  weather  prediction, 
was  tried.  This  first  system  effectively  took  a gradient  and  multiplied  it  by 
a speed,  although  what  was  really  programmed  was  an  interpolation  of  the  energy 
a distance  upstream  depending  on  the  component  group  velocity.  This  does  not 
work  because  in  each  future  time  step  the  energy  spreads  further  and  further 
as  illustrated  schematically  by  the  solid  vertical  lines  in  Figure  3-2.  Note, 
for  instance,  that  these  vertical  lines  decrease  in  amplitude  during  each  time 
step,  and  new  lines  form  further  and  further  away  from  the  hatched  area,  which 
is  where  they  should  be.  Higher  order  interpolation  does  not  help.  The  cor- 
rect method  must  allow  for  discontinuities.  These  sharp  changes  are,  of  course, 
often  observed  with  the  passage  of  a cold  front  and  with  the  onset  of  swell. 

This  is  a major  difference  between  the  present  effort  and  that  of  Gelci  and 
Chavy  (Reference  16,  1961). 

A simple  approach  was  used  to  establish  a method  of  propagation,  though  others 
can  be  devised  with  more  accuracy.  In  principle,  what  was  done  was  to  keep 
track  of  how  far  past  and  to  the  side  of  an  average  grid  point  the  wave  com- 
ponents were,  then  when  the  component  had  moved  far  enough  to  reach  the  ad- 
jacent grid  point,  the  energy  was  "jumped".  Since  the  grid  system,  in  use,  is 
only  approximate,  this  was  simplified  by  keeping  only  one  record  for  each  com- 
ponent for  all  of  the  grid  points,  because  at  all  grid  points  the  same  compo- 
nent moves  in  the  same  direction  at  the  same  speed.  Thus,  if  it  takes,  say, 
three  time  steps  for  a particular  component  to  travel  to  the  adjacent  grid 
point  on  every  third  step,  the  entire  field  will  be  jumped.  That  is  to  say 
the  energy  from  grid  point  upstream  replaces  the  energy  at  the  adjacent  grid 
point  downstream. 

Because  of  the  grid  shape,  north-south  and  east-west  components  are  very  simple. 
However,  in  other  directions  it  is  necessary  to  remember  how  far  the  wave  com- 
ponent exists  to  the  sides  as  well  as  forward  and  backward.  Here,  the  choice 
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ASSUME  (I)  PROPAGATION  VELOCITY  EQUALS  1/2 
THE  DISTANCE  BETWEEN  GRID  POINTS 
(2)  NO  CHANGE  IN  ENERGY  OF  THE  COMPONENTS 


Figure  3-2  Results  of  Propagating  Wave  Energy,  Schematic 
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of  30  degrees  directional  increments  was  fortuitous.  If  an  energy  pocket 
travels  the  maximum  distance  within  the  chosen  grid,  it  will  only  he  off  hy 
one  grid  point  if  it  is  Jumped  alternately  toward,  say,  West,  then  Southwest. 
What  is  actually  programmed  is  that  a Jump  is  made  in  the  cardinal  or  inter- 
mediate direction,  depending  on  whether  the  time  step  number  is  even  or  odd 
when  the  energy  packet  has  traveled  an  average  distance.  The  method  of  pro- 
pagation used  can  have  a strong  effect  on  the  shape  and  size  of  the  grid.  With 
four  individual  distance  markers  for  each  frequency-directional  component,  it 
might  be  possible  to  increase  the  grid  size  and  increase  the  accuracy  of  the 
computations  at  the  same  time. 


3-5  EQUATIONS 

The  Neumann  non-directional  spectrum,  along  with  the  SWOP  directional  effects, 
were  taken  as  the  maximum  spectrum  for  a given  wind  condition.  This  is  de- 
fined as  the  finite  difference  equation: 
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AS  = maximum  energy  in  the  increment, 
max 


41  *■ 


The  growth  function  from  the  CCS  curves  is  also  needed.  No  simple  curve 
could  be  fitted,  so  this  was  tabulated  as  presented  in  Table  3-1  and  used 
with  linear  interpolation  for  intermediate  values.  These  values  were  read 
directly  from  the  CCS  curves  presented  in  H.O.  603  for  two  hour  increments, 
then  smoothed  slightly.  Note  that  the  increase  in  energy  in  each  two  hour 
time  increment  depends  on  both  the  wind  speed  and  the  initial  energy  state. 

The  definition  of  this  initial  energy  state  was  made  arbitrarily.  To  be  con- 
sistent with  the  CCS  curves  for  constant  wind  conditions,  it  was  decided  to 
find  this  initial  total  energy  by  adding  all  components  within  ±90  degrees 
from  the  direction  of  the  wind.  The  consequences  of  this  assumption  are  that 
long  swell  is  added  in  the  same  way  as  local  sea.  This  means  that  when  swell 
is  present  with  a similar  direction,  the  sea  will  build  up  faster  than  other- 
wise. Another  effect  is  that  the  more  changeable  the  wind  direction,  the 
more  slowly  the  waves  will  grow. 

Instead  of  these  assumptions,  the  second  parameter  of  the  growth  function 
table  could  have  been  a mean  frequency  or  an  upper  frequency.  In  this  case, 
some  more  complicated  assumptions  would  have  been  required  regarding  swell 
and  wind  direction  changes. 


3-6  COMPUTATIONAL  STEPS 

With  all  of  the  above  methods,  assumptions,  and  equations,  it  is  possible  to 
present  the  major  individual  steps  of  the  computation.  A master  flow  diagram 
and  program  listing  are  presented  in  Appendices  A and  B,  respectively.  These 
steps  are  repeated  for  each  time  step,  as  follows: 

1.  Given  initial  spectrum  specified  by  the  120  components  at  each  grid 
point . 
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Table  3-1  Energy  Growth  in  Feet  of  the  Wave  Spectrum 
Per  Two  Hour  Time  Step 


(ft2) 

10 

14 

18 

22 

26 

30 

34 

38 

44 

50 

56 

0 

0.1 

0.15 

0.2 

0.35 

0.5 

0.7 

1.0 

1-5 

2.0 

3-1 

4.5 

1 

0.0 

0.6 

1.7 

1.8 

2.0 

2.5 

3-0 

3-5 

4.0 

5-3 

7-0 

5 

0.0 

0.0 

3.2 

3-4 

3-8 

4.0 

4.5 

6.0 

8.7 

12.0 

10 

2.0 

4.4 

4.5 

4.5 

5.0 

7-5 

11.0 

15.0 

15 

0.0 

5-0 

5-0 

5-0 

5.4 

8.0 

12.0 

17-0 

20 

5.0 

5-5 

6.0 

6.4 

9-0 

13-0 

20.0 

25 

4.5 

5-9 

6.4 

6.9 

9-4 

14.0 

23-0 

30 

0.0 

6.6 

7-0 

8.0 

10.0 

15-0 

25-0 

35 

7-7 

8.0 

9-0 

10.4 

16.0 

25.0 

40 

8.7 

9-1 

9-8 

10.7 

17.0 

25-0 

45 

9-0 

10.0 

10.5 

11.0 

17-5 

25.0 

50 

7.0 

11.0 

11.0 

11.4 

18.0 

25.0 

60 

0.0 

12.0 

12.0 

12.0 

18-5 

25.0 

70 

12.0 

12.4 

12.4 

19.0 

25.0 

80 

12.0 

12.6 

12.8 

20.0 

25.0 

90 

11.0 

12-9 

13.2 

21.0 

25-0 

100 

10.0 

13-2 

13.8 

22-0 

25.0 

115 

0.0 

13.5 

14.2 

23-0 

25.0 

130 

14.0 

15.0 

24.0 

25.0 

145 

14.0 

15.8 

25.0 

25.0 

160 

13-0 

17.0 

25.0 

25.0 

175 

10.0 

20.0 

25-0 

25-0 

200 

0.0 

23.0 

25-0 

25.0 

/v 
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2.  Given  wind  velocity  at  each  grid  point  for  the  succeeding  time  interval. 

3.  Sum  the  total  spectral  energy  at  each  grid  point  within  ±90  degrees  of 
the  wind  direction. 

4.  Look  up  the  total  allowable  change  in  energy  for  each  grid  point  from  the 
growth  function  tabulation. 

5.  Find  the  maximum  possible  spectrum  at  each  grid  point  from  equation  (l). 

6.  Compare  the  initial  spectrum  with  the  maximum  spectrum  for  each  grid 
point  and  note  any  gaps  that  can  be  filled  in. 

7-  Spread  the  available  growth  energy  into  the  available  gaps  starting  with 
the  high  frequencies  first  for  each  grid  point.  Secondarily,  the  direc- 
tions nearest  to  that  of  the  wind  are  to  be  used  first. 

NOTE:  This  procedure  finds  the  total  energy 
growth  to  be  expected,  then  splits  it  into 
component  directions  and  frequencies.  The 
method  carefully  includes  any  higher  frequencies 
which  may  have  been  missing  in  the  initial 
spectrum. 

8.  The  final  step  is  to  allow  each  component  of  the  spectrum  to  propagate  in 
its  own  direction  and  at  its  own  group  velocity.  The  spectrum  at  each 
grid  point  is  then  the  initial  condition  for  the  next  time  step. 


3-7  INITIAL  SPECIFICATION 

In  order  to  specify  the  spectrum  for  the  initial  conditions  properly  in  the 
first  time  step,  120  x 519  = 62,280  4E  values  should  be  read  into  the  com- 
puter. If  all  these  numbers  were  available  on  maps,  the  task  of  scaling 
them  off  would  be  almost  impossible.  However,  present  spectral  observations 
are  so  sparse  that  the  initial  spectral  fields  are  not  available.  Instead, 
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a very  rough  first  guess  is  used.  Estimates  of  the  significant  height  and 
earlier  winds  are  used  to  get  a rough  estimate  of  the  initial  spectrum.  This 
is  accomplished  by  assuming  for  each  grid  point  that  the  previous  wind 
(at  t = -l)  had  been  constant  long  enough  to  raise  the  specified  sea.  If 
the  wind  was  not  high  enough  to  raise  the  sea,  the  excess  spectral  energy 
was  discarded. 

There  are  so  many  errors  in  this  process  that  they  are  not  worth  detailing. 
However,  after  many  time  steps,  all  of  this  initial  energy  will  be  dissipated 
or  masked  by  the  effects  of  later  winds  making  the  remaining  errors  insigni- 
ficant. Thus,  the  first  few  time  steps  cannot  give  correct  answers.  This 
means  that  any  forecast  must  begin  with  a hindcast  to  find  the  present  spec- 
tral conditions.  Present  limited  experience  suggests  about  one  day  for  this 
hindcast. 

3-8  ILLUSTRATIVE  COMPUTATION 

In  order  to  see  how  the  model  behaves  with  this  computing  scheme  an  artificial 
test  case  was  run.  The  initial  significant  height  to  start  the  computation 
was  assumed  to  be  20  feet.  Then,  a constant  wind  of  47  knots  from  270  degrees 
was  allowed  to  blow  over  only  this  same  grid  point  for  9 time  steps  (18  hours) 
as  well  as  for  the  initial  step.  The  results  of  this  experiment  are  shown  in 
Figure  3-3  along  with  a tenth  time  step  having  no  wind.  For  comparison,  a 
47  knot  wind  can  generate  waves  having  a significant  height  of  67  feet  for 
infinite  fetch  and  duration.  For  fetches  of  120  nm,  the  significant  height 
should  be  about  21  feet.  It  should  also  be  pointed  out  that  because  of  non- 
linear growth  characteristics,  the  final  forecast  height  is  a function  of  the 
initial  conditions. 

This  impossible  case  illustrates  the  possible  accuracy  of  the  model  by  show- 
ing the  "minor"  variations  resulting  from  the  discrete  and  discontinuous 
assumptions . When  energy  is  propagated  by  jumps,  the  energy  at  the  location 
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it  has  left  is  decreased  and  vice-versa.  If  the  propagation  phase  were  to 
be  carried  out  at  the  beginning  of  each  time  step  instead  of  as  the  final 
phase,  these  variations  would  be  more  masked  when  the  seas  were  near  full 
development.  But  they  would  still  affect  the  accuracy  of  succeeding  steps 
in  the  same  way. 

This  example  shows  what  the  system  would  predict  near  the  rear  boundary  of  a 
classical  sharp  edged  fetch.  However,  in  the  real  case,  there  are  usually 
some  smaller  waves  behind  the  fetch.  The  energy  from  these  smaller  waves 
also  propagates  into  the  fetch  to  replace  the  energy  propagating  out  and 
decreasing  the  irregularities  illustrated  here.  Also,  any  grid  point  within 
the  fetch  would  receive  energy  from  the  surrounding  grid  points  making  these 
errors  relatively  insignificant. 


Several  other  points  are  illustrated  in  this  example.  Note  that  the  spectrum 
is  not  completely  symmetrical  because  odd  amounts  of  energy  have  not  been 
split  symmetrically  but,  rather,  used  in  a particular  order  so  that  a slight 
bias  toward  one  side  exists.  If  later  study  shows  this  to  be  important,  the 
program  can  be  corrected  easily.  Another  point  of  interest  is  the  magnitude 
of  the  swell  increments  that  move  out.  They  are  governed  by  the  amount  of 
spectral  energy  in  each  frequency-directional  component . Thus,  since  spectra 
have  their  peaks  near  the  low  frequency  cut  off,  smaller  frequency  increments 
must  be  used  in  the  lower  frequencies  to  keep  the  jumps  from  being  too  large. 
This  was  the  reason  the  frequency  intervals  were  chosen  to  have  midpoints 
spaced  at  constant  period  increments . Experience  may  well  show  that  even 
more  drastic  biasing  toward  low  frequencies  is  better.  Lastly,  it  should  be 
pointed  out  that  the  swell  shown  is  a function  of  the  time  step  number  on 
which  the  computation  was  begun.  By  starting  at  time  step  zero,  no  jumps 
took  place  until  time  step  4.  Starting  at  some  other  step  could  have  caused 
the  jumps  to  take  place  sooner  and  in  a different  order.  This  is  because  the 
distance  that  the  energy  of  the  component  has  traveled  past  a grid  point  is 
treated  as  a function  of  the  time  step  count  only. 


3-14 

LOCKHEED  MISSILES  & SPACE  COMPANY 


IMSC  -801296 


CHAPTER  4 

HINDCASTS  AND  VERIFICATION 

4-1  DISCUSSION 

To  maintain  objectivity  and  mechanize  preparation  of  the  required  input  data 
for  the  IBM  7090,  rules  shown  in  Table  4-1  were  established. 


Table  4-1  General  Instructions  for  Preparing  Input  Data 


Item 


Description 


6 

7 


All  surface  wind  speeds  used  are  averages  over  the  square  centered 
at  the  grid  point  with  sides  of  120  nautical  miles. 

At  time  t=0,  the  significant  wave  height  at  every  grid  point  is 
needed.  Directions  are  not  used. 

Surface  winds  are  needed  from  the  6-hourly  map  previous  to  t=0  up  to 
at  least  t=48  hours  in  6-hour  or  shorter  steps.  If  information  is 
available  at  less  than  6-hourly  steps,  it  can  be  used  down  to  2 hours. 

Wind  speeds  less  than  10  kts  can  be  recorded  as  0.  All  omissions  are 
considered  the  same  as  0. 

Direction  is  to  be  recorded  in  relation  to  the  grid,  not  true  North. 
The  top  of  the  grid  is  North.  These  two  values  will  be  the  same  at 
35 # longitude. 

Direction  from  the  "North"  can  be  recorded  as  either  "000"  or  "360". 

The  speed  is  recorded  first,  then  the  direction  in  3 digits, 

(as  27.020.). 
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Figure  4-1  gives  the  location  number  assigned  to  each  of  the  grldpolnts. 

These  numbers  are  used  for  specifying  the  location  where  a complete  spectrum 
is  desired  and  for  many  operations  internal  to  the  machine.  The  IBM  7090 
program  required  that  the  wind  speeds  or  wave  heights  be  punched  Into  stand- 
ard Hollerith  cards  with  18  four-place  values  per  punch  card.  For  winds,  two 
four-place  numbers  are  needed  to  specify  each  grid  point  having  two  values 
each  for  only  nine  grid  points  on  each  card.  Thus,  58  cards  are  required  to 
specify  the  winds  at  each  two-hour  time  step. 

When  this  study  was  begun,  no  adequate  ocean  wave  spectra,  from  measured  waves 
in  severe  sea  conditions,  were  available.  Therefore,  use  of  synoptic  map6  of 
significant  height  was  planned.  In  order  to  be  as  objective  as  possible,  the 
U.  S.  Navy  Hydrographic  Office  was  asked  to  choose  a severe  storm  condition 
and  to  provide  copies  of  their  regular  synoptic  wave  charts.  Careful  stream- 
line analyses  were  made  of  the  ten  six-hourly  synoptic  charts,  but  the  wave 
charts  were  accepted  without  change  to  maintain  objectivity.  The  period  chosen 
was  8 through  10  September  1956.  Later,  during  the  course  of  this  research, 
waves  were  measured  for  a severe  sea  condition  during  l6  December  through 
18  December  1959 , from  which  spectra  were  computed.  This  case  is  presented 
here  as  a second,  and  perhaps  more  meaningful,  test  of  the  model  established. 
The  synoptic  situations,  input  data,  results  and  verifications  of  the  two 
cases  are  presented  and  discussed  in  a succeeding  separate  section. 


4-2  CASE  I,  8 THROUGH  10  SEPTEMBER  1956 

There  were  three  important  aspects  of  the  synoptic  situation  for  the  North 
Atlantic  Ocean  at  the  beginning  of  this  period.  First,  a cold  front  was 
roughly  parallel  to  the  East  Coast  of  North  America  and  had  moved  200  to  300 
miles  into  the  ocean.  This  cold  front  extended  from  an  occlusion  and  deep 
cyclone  near  the  southern  part  of  Greenland.  On  the  southern  extremity,  the 
front  was  relatively  weak  and  tending  to  become  stationary  with  the  remains 
of  a hurricane  approximately  4o0  miles  ahead  of  the  front.  Another  deep  low 
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Figure  4-1  Number  Location  Assigned  to  the  Various  Grid  Points 
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was  sitting  about  5 00  miles  off  the  coast  of  France,  while  the  central  part  of 
the  ocean  was  dominated  by  the  usual  anticyclone.  The  analysis  made  by  the 
National  Weather  Analysis  Center,  is  shown  in  Figure  4-2(a).  The  Initial  sea- 
condition  chart  prepared  by  the  U.  S.  Navy  Fleet  Weather  Central  is  presented 
in  Figure  4-3(a) . The  remaining  parts  of  Figures  4 -2(a)  through  4-2(f ) and 
Figures  4-3(a)  through  4-3(c)  show  analyses  for  the  succeeding  times.  During 
these  succeeding  periods,  the  cyclone  off  France  deepened  and  remained  fairly 
stationary,  providing  northerly  winds  of  about  40  knots  over  what  might  be 
considered  a 300  mile  fetch  and  several  days  duration.  In  turn,  the  weak  cy- 
clone off  Florida  became  a wave  on  the  front  and  deepened  to  provide  a compli- 
cated moving  fetch  with  wind  speeds  up  to  50  knots. 

The  initial  significant  height  of  the  waves  at  each  grid  point  was  interpolated 
from  the  initial  sea  condition  chart.  Similarly,  the  synoptic  charts  were  used 
to  prepare  surface  streamline  charts  that  were  then  scaled  in  accordance  with 
the  rules  of  Table  4-1.  The  data  actually  used  is  presented  in  Appendix  C. 
Hindcasts  were  prepared  from  these  winds  and  initial  conditions  and  the  signifi- 
cant height  field  for  every  time  step  was  tabulated  along  with  many  individual 
spectra.  An  average  of  about  two  minutes  per  time  step  was  required  for  the 
machine  computation.  However,  only  the  initial  field  and  results  at  12-hour 
intervals  are  presented  in  Figure  4 -4 (a). 

Comparison  of  the  significant  height  fields  in  Figures  4-3(a)  through  4-3(c), 
and  Figures  4-4(a)  through  4-4(e)  show  fairly  good  agreement  except  in  the  high 
wave  area  near  the  European  coast.  However,  with  respect  to  this  area,  a 
graphical  forecast  using  the  methods  of  H.O.  603  would  have  given  values  close 
to  those  found  by  the  machine  forecast.  This  and  several  ship  wave  reports 
confirmed  the  high  wave  area.  Therefore,  an  independent  verification  was 
attempted  using  the  estimated  wave  heights  recorded  by  the  ocean  weather  ships. 
In  order  to  check  the  wave  charts  and  get  a large  enough  sample  to  provide  a 
reasonable  comparison,  information  from  all  available  weather  ships  from  the 
files  of  the  National  Weather  Records  Center  of  the  U.  S.  Weather  Bureau  were 
used.  Plots  of  the  comparison  for  each  ship  are  shown  in  Figure  4-5.  Ship  J, 
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which  was  in  the  high  wave  region  off  Europe,  confirmed  the  inadequacy  of  the 
operational  wave  charts  for  verification  purposes.  No  comparison  of  the  periods 
was  made  because  of  the  well  known  inadequacy  in  the  methods  of  observing  and 
recording  periods. 

These  verification  graphs  in  Figure  4-5  show  that  the  model  forecast  the 
heights  about  as  well  as  the  operational  wave  charts.  The  major  failing  of 
the  model  is  the  continued  forecasting  of  high  seas  after  the  seas  have  lowered. 
This  type  of  error  also  occurred  in  the  Case  II  forecast  and  is  discussed  in  de- 
tail in  a separate  section. 

In  this  test,  the  method  seems  quite  good  except  for  this  lack  of  dissipation. 
The  important  result  is  that  large  scale  patterns  can  be  represented  easily. 


4-3  CASE  II,  16  THROUGH  18  DECEMBER  1959 

All  the  data  for  both  the  wind  fields  and  the  waves  for  this  case  were  provided 
through  the  courtesy  of  New  York  University,*  as  shown  in  Appendix  4.  This 
particular  storm  is  being  studied  by  a number  of  other  scientists,  as  described 
in  Appendix  5-  The  synoptic  situation  and  general  forecast  problem  is  described 
by  L.  Moskowitz  as  follows;  (W.J.  Pierson,  Jr.,  1961,  Personal  Communication) 

"The  week  beginning  on  15  December  1959  saw  several  frontal  systems  pass  thr  gh 
that  area  of  the  North  Atlantic  with  which  this  paper  deals.  Position  J,  sit- 
uated at  approximately  52o40'N  and  20°W,  lies  fax  to  the  north  of  the  subtropi- 
cal anticyclone  and  lies  south  of  the  major  tracks  of  the  North  Atlantic  cyclones. 
Storms  in  this  area  of  the  North  Atlantic  are  known  for  their  persistence 
and  intensity.  The  most  impressive  aspect  of  these  storms  is  the  development 
of  the  state  of  the  sea. 


*See  paragraph  2 of  the  Foreword. 
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"From  the  point  of  view  of  the  wave  forecaster  the  storm  of  l6  - l8  December  is 
of  interest.  The  storm  developed  from  a wave  cyclone  which  had  formed  on  a cold 
front  on  15  December  at  about  12Z.  Between  18Z  on  the  15th  and  00Z  on  the  l6th, 
the  storm  had  developed  its  own  circulation  and  the  central  pressure  dropped  to 
984mb  at  00Z  on  the  16th.  At  this  time  the  frontal  system  associated  with  the 
storm  had  partially  occluded.  The  storm  center  moved  rapidly  towards  the  north- 
east and  on  the  16th,  at  about  l8z,  had  a central  pressure  of  960mb  and  was  lo- 
cated at  57-5°N  and  21°W.  Winds  behind  the  storm  on  the  order  of  fifty  knots 
had  developed  between  06  and  12Z  on  the  l6th  and  persisted  through  l8z  on 
17  December.  At  about  00Z  on  the  17th,  the  storm  center  began  to  stagnate  at 
about  60°N  and  20°W,  however,  the  frontal  system  continued  its  rapid  eastward 
movement.  At  about  l8z  on  the  17th  the  storm  center  had  begun  to  fill.  At 
06Z  on  the  18th  the  central  pressure  had  already  increased  to  972mb. 


"At  18Z  on  l6  December,  the  OWS  Weather  Reporter  passed  through  the  occluded 
front  associated  with  the  storm  center.  At  this  time  thirty  knot  winds  were 
reported.  Six  hours  later  the  winds  had  built  up  to  60  knots  from  the  south- 
west. The  significant  wave  height,  jumped  from  15-7'  at  l8z  on 

16  December  to  26.2'  at  00Z  on  the  17th.  Winds  on  the  order  of  sixty  knots 
(Beaufort  ll)  were  not  uncommon  in  the  general  area  during  the  above  times. 

It  is  to  be  noted  that  all  wind  data  until  12Z  on  17  December  were  Beaufort 
estimates.  The  maximum  of  39-7'  occurred  at  18Z  on  17  December.  The 

wind  speed  at  this  time  was  48  knots.  The  wave  heights  began  to  decrease 
steadily  thereafter,  except  for  a secondary  maximum  of  35*4'  at  00Z  on 
18  December. 


"When  the  data  station  (OWS  in  this  case)  is  not  stationary,  the  forecasting 
problem  is  further  complicated.  Such  is  this  case.  At  00Z  on  15  December,  the 
OWS  at  position  J was  already  underway  for  its  home  port.  The  OWS  "Weather 
Reporter",  the  relief  ship,  was  also  underway  heading  for  position  J at  the 
same  time.  Tne  two  ships  passed  each  other  sometime  between  00Z  and  06Z  on 
16  December  (presumed  from  synoptic  charts).  Since  the  wave  records  used  were 
taken  from  the  OWS  "Weather  Reporter",  the  fetches  to  be  used  in  forecasts 
Decome  even  more  critical.  Not  until  l8Z  on  17  December  did  the  "Weather 
Reporter"  reach  position  J even  though  she  reported  she  was  at  position  J as 
early  as  12Z  on  17  December." 
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Hindcasts  were  made  for  the  entire  area,  but  only  the  significant  heights  for 
the  area  near  the  observing  ship  are  shown  in  Figures  4-6(a)  through  4- 6(f). 

The  crosses  on  this  Figure  show  the  relative  position  of  the  observing  ship  at 
the  respective  times.  The  complete  spectra  for  nearby  grid  points  was  also 
printed.  Some  examples  of  this  are  shown  in  Tables  4-2 (a)  through  4-2 (c). 
Because  of  the  slight  difference  in  the  reporting  and  computing  intervals, 
a comparison  chart  for  verification  of  significant  heights  was  plotted  and 
is  presented  as  Figure  4-7. 

"’here  are  three  significant  points  shown  in  Figure  4-7  that  require  discussion. 

1 >st  of  the  11-17  foot  significant  heights  in  the  initial  observations  were 
due  to  swell,  because  the  winds  were  only  11  - 17  knots.  The  forecasting  model 
drops  initial  swell  so  that  it  automatically  cut  this  to  about  two  feet. 

Despite  this  expected  initial  handicap  and  the  rapidly  rising  sea,  the  fore- 
casting model  caught  up  with  the  observed  35  foot  sea  after  only  one  day  of 
hindcasting,  and  followed  the  observed  significant  height  in  the  succeeding 
twelve  hours.  As  seen  in  Figure  4-7,  the  forecast  rise  in  height  is  parallel 
to  that  observed  with  a time  lag  due  to  poor  initial  conditions . At  this 
time,  the  winds  dropped  and  changed  direction  and  the  sea  was  observed  to 
lower  rapidly  to  below  15  feet  in  twelve  hours . The  model  did  not  predict 
this  rapid  decrease.  Possible  reasons  for  this  error  are  discussed  in  the 
following  section. 

A comparison  of  the  observed  and  forecast  values  of  the  frequency  spectrum  was 
made  with  disappointing,  but  not  unsuspected  results,  as  shown  in  Figures  4-8(a) 
through  4-8(c).  In  these  figures,  the  plotted  points  represent  the  estimated 
spectra  and  the  90$  confidence  interval  integrated  over  a frequency  increment  of 
0.00555  sec”1.  The  stippled  regions  show  the  forecasts  on  the  same  scale.  The 
model  forecasts  do  not  have  as  much  low  frequency  energy  as  the  spectra  computed 
from  the  observations  of  the  OWS  Weather  Reporter.  There  are  several  possible 
causes  for  this  discrepancy.  They  are  discussed  in  the  following  section. 
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Figure  J*-6(a)  Forecast  Maps  at  the  End  of  the  Respective  Time  Steps 
and.  the  Observed  Significant  Heights  for  Verification 
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Figure  4 -6(b)  Forecast  Maps  at  the  End  of  the  Respective  Time  Steps 
and  the  Observed  Significant  Heights  for  Verification 
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Figure  4-6(c)  Forecast  Maps  at  the  End  of  the  Respective  Time  Steps 
and  the  Observed  Significant  Heights  for  Verification 
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Figure  4-6(d)  Forecast  Map6  at  the  End  of  the  Respective  Time  Steps 
and  the  Observed  Significant  Heights  for  Verification 
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Figure  U-6(e)  Forecast  Maps  at  the  End  of  the  Respective  Time  Steps 
and  the  Observed  Significant  Heights  for  Verification 
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Table  4-2(a)  Example  Forecast  Directional  Spectrum 


SPECIAL 
CASE  2 

TIME  = STEP  12 . GRIDPOINT  = lb 


9=0 

0.2203E-01 

0.4306E-00 

0.1447E  01 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

9 = 30 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

9 = 60 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

9 = 90 

0.1196E-01 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

9 = 120 

0.3515E-01 

0.7476E-01 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

9 = 150 

0.5865E-01 

0.3076E-00 

0.2273E  01 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

9 = l80 

0.7055E-01 

0.1543E-01 

0.5062E  01 

0.8767E  01 

0. 

0. 

0. 

0. 

0. 

0. 

9 = 210 

0.9283E-01 

0.1702E  01 

0.5493E  01 

0.6646E  01 

0.8355E  00 

0. 

0. 

0. 

0. 

0. 

O 

CVl 

II 

CD 

0.1235E-00 

0.2297E  01 

0.6446E  01 

0.1235E  02 

0.1562E-01 

0.9765E-03 

0. 

0. 

0. 

0. 

9 = 270 

0.1291E-00 

0.2515E  01 

0.8537E  01 

0.1826E  02 

0.8789E-02 

0. 9765E-03 

0. 

0.2691E-01 

0.2990E-02 

0. 

9 = 300 

0.1076E-00 

0.2088E  01 

0.6727E  01 

0.1235E  02 

0. 

0. 

0. 

0. 

0. 

0. 

9 = 330 

0.6250E-01 

0.1248E  01 

0.4189E  01 

0.1318E-01 

0. 

0. 

0. 

0. 

0. 

0. 

For  each  direction,  frequency  coordinates  are  as  follows: 

0.35  0.188  0.129  O.O98  0.079 

0.066  0.057  0.050  0.044  0.04 


—1  2 

Units  are  (deg.),  (sec-  ),  and  (ft  /frequency  interval),  respectively. 


Spectral  energies  are  presented  in  standard  machine  language 
"floating  point"  format,  in  which  the  numbers  and  signs 
following  the  Es  are  the  exponents  of  the  base  10  by  which 
the  characteristics  must  be  multiplied. 
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Table  4-2 (b)  Example  Forecast  Directional  Spectrum 


SPECIAL 
CASE  2 

TIME  = STEP  15.GRIDP0INT  = 6l 


e = o 

0.7910E-01 

0.1621E  01 

0.5275E  01 

0.5420E  01 

0. 

0. 

0. 

0. 

0. 

0. 

e = 30 

0.2203E-01 

0.6857E  00 

0.1196E  01 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

S 

II 

CD 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

© = 90 

0.1977E-01 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

© = 120 

0.1977E-01 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

e = 150 

O.5865E-OI 

O.3O76E-OO 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

© = 180 

0.7910E-01 

0.8643E  00 

0.1999E  01 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

© = 210 

0.1235E-00 

0.1543E  01 

0.4189E  01 

0.1128E  01 

0.6103E-04 

0. 

0. 

0. 

0. 

0. 

0 = 240 

0.1235E-00 

0.2226E  01 

0.7097E  01 

O.4785E  01 

0.1562E-01 

0.2197E-02 

0. 

0.1406E-00 

0.2441E-01 

0. 

<D 

II 

ro 

-v] 

0 

0.1291E-00 

0.2640E  01 

0.6606s  01 

0.1092E  02 

0. 

0.1525E-02 

0. 

0. 

0. 

0. 

© = 300 

0.1291E-00 

0.2615E  01 

0.8491E  01 

0.1066E  02 

0.1983E-02 

O.I763E-OI 

0. 

0.2441E-03 

0. 

0. 

© = 330 

0.1235E-00 

0.2465E  01 

0.7998E  01 

0.1714E  02 

0. 

0. 

0. 

0. 

0. 

0. 

For  each  direction,  frequency  coordinates  are  as  follows: 

0.35  0.188  0.129  0.098  0.079 

0.066  0.057  0.050  0.044  0.04 


Units  are  (deg.),  sec"1),  and  ( ft*/ frequency  interval),  respectively. 


Spectral  energies  are  presented  in  standard  machine  language 
"floating  point"  format,  in  which  the  numbers  and  signs 
following  the  Es  are  the  exponents  of  the  base  10  by  which 
the  characteristics  must  be  multiplied. 


4-30 


VJ 


I 


LOCKHEED  MISSILES  ft  SPACE  COMPANY 


CHANGE  1 

imsc -801296 


Table  4— 2(c)  Example  Forecast  Directional  Spectrum 


SPECIAL 
CASE  2 

TIME  = STEP  18.GRIDP0INT  = 73 


9 = 0 

0.8355E-01 

0.1621E-01 

0.5275E  01 

0.1112E  02 

0.7097E  01 

0. 

0. 

0. 

0. 

0. 

9 = 30 

0.4785E-01 

0.5742E  00 

0.2615E  01 

0.2226E  01 

0. 

0. 

0. 

0. 

0. 

0. 

9 = 60 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

9 = 90 

0.1196E-01 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

9 = 120 

0.1562E-01 

0.7476E-01 

0.7656E  00 

0.2658E-00 

0. 

0. 

0. 

0. 

0. 

0. 

9 = 150 

0.4785E-01 

0.1283E  01 

0.2197E-00 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

<D 

II 

H 

CO 

O 

0.4785E-01 

0.0999E  01 

0.1485E  01 

0.7793E  00 

0. 

0. 

0. 

0. 

0. 

0. 

9 = 210 

0.9765E-01 

0.1466E  01 

0.4189E  01 

0.2226E  01 

0.1428E  01 

0. 

0. 

0. 

0. 

0. 

9 = 240 

0.1235E-00 

0.2088E  01 

0.7097E  01 

0.1501E  02 

0.9474E  01 

0.3076E-00 

0.5493E-03 

0. 

0. 

0. 

O 

CM 

II 

a> 

0.1348E-00 

O.2615E  01 

O.813IE  01 

0.1760E  02 

0.1812E  02 

0. 

0.1525E-02 

0. 

0. 

0. 

9 = 300 

0.1291E-00 

0.2640E  01 

O.8767E  01 

0.1575E  02 

0.3315E  02 

0.8767E  01 

0. 

0. 

0. 

0. 

<D 

II 

UJ 

CO 

0 

0.1181E-00 

0.2392E  01 

0.7998E  01 

0.1235E  02 

0.1839E  02 

0. 

0. 

0. 

0. 

0. 

For  each  direction,  frequency  coordinates  are  as  follows: 

0.35  0.188  0.129  0.098  0.079 

0.066  0.057  0.050  0.044  0.04 

2 

Units  are  (deg.),  sec"  ),  and  (ft  /frequency  interval),  respectively. 

Spectral  energies  are  presented  in  standard  machine  language 
"floating  point"  format,  in  which  the  numbers  and  signs 
following  the  Es  are  the  exponents  of  the  base  10  by  which 
the  characteristics  must  be  multiplied. 
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DECEMBER  16  DECEMBER  17  DECEMBER  18 

- 18  - - 00  - - 06  - - 12  - - 18  - - 00  - - 06  - 

ABSCISSA  SHOWS  STEP  NUMBER,  VERIFICATION  WAVE  OBSERVATION 
USEu,  AND  THE  WIND  OBSERVATION  TIME  USED  IN  THE  STEP 


Figure  U~7  Case  II,  Verification  Graphs 
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CHAPTER  5 
CONCLUSIONS 


/ 


C 


5-1  SUMMARX  OF  RESULTS 

The  -results  described  in  the  preceding  section  are  promising  and  would  be 
excellent  except  for  three  conditions.  They  are:  l)  too  few  low  frequency 
waves  in  the  spectrum,  2)  two  slow  a decay  for  swell,  and,  3)  the  lack  of 
knowledge  about  the  initial  conditions. 

It  has  been  shown  that,  l)  numerical  forecasts  of  sea  and  swell  spectra  for 
the  entire  ocean  can  be  made,  2)  wave  forecasts  can  be  made  without  defining 
a fetch,  thus  decreasing  the  subjectivity  of  the  forecast,  and  3)  a completely 
objective  model  has  been  established  which  should  be  useful  for  testing  hypoth- 
eses and  upon  which  further  research  can  be  based.  This  model  keeps  track  of 
all  frequencies  involved  and  propagates  them  with  the  correct  group  velocities. 


5-2  DISCUSSION  OF  ERRORS 

There  are  several  possible  causes  for  having  too  few  low  frequency  waves  in  the 
spectrum,  under  the  assumption  that  the  observations  are  correct.  Long  swell, 
which  the  model  would  ignore,  could  have  been  present  at  some  distance  away  in 
the  initial  sea.  However,  this  is  unlikely  because  of  the  smooth  shapes  of 
the  observed  spectra.  It  was  assumed  that  the  energy  increment  added  was  a 
function  of  the  energy  within  ±90  degrees  of  the  wind  direction.  If  this  is 
changed  to  consider  energy  in  all  directions,  lower  frequencies  and  higher 
waves  will  result.  A more  likely  possibility  is  simply  that  the  Neumann  spec- 
trum does  not  predict  energy  at  low  enough  frequency.  This  has  been  pointed 
out  in  another  study.  (W.  J.  Pierson,  Jr.,  Reference  12) 
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Another  Important  consideration  is  the  fact  that  the  Neumann  spectrum  was  not 
devised  for  this  particular  purpose.  It  assumes  a constant  wind  throughout 
the  length  of  the  fetch.  However,  in  the  real  case  which  the  present  model 
allows,  the  wind  direction  often  has  considerable  variation  in  time.  Such  a 
changeable  wind  direction  tends  to  spread  the  higher  frequencies  into  other 
directions  before  allowing  the  lower  frequencies  to  grow.  A reasonable  correc- 
tion for  this  variation  might  be  a small  energy  transfer  function  which  acted 
to  transfer  energy  from  that  established  by  previous  wind  directions  to  that 
which  is  consistent  with  the  latest  wind  velocity. 

Swell  did  not  decay  fast  enough  in  either  of  the  cases  studied.  This  can  be 
caused  by  either  the  lack  of  a dissipation  function  or  by  improper  dispersion. 
Because  of  the  discrete  definition  of  the  spectrum,  with  all  of  the  energy  with- 
in a 30  degree  sector  concentrated  along  the  central  direction,  waves  from  a 
single  grid  point  are  not  allowed  to  disperse  properly.  When  the  swell  compo- 
nent has  traveled  a distance  of  about  three  grid  points  from  the  location  of 
its  generation  the  value  of  the  component  should  have  decreased  by  approxi- 
mately one -half.  The  remaining  half  of  the  energy  should  have  spread  into 
the  surrounding  grid  points  which  are  in  the  triangle  between  the  30  degree 
radial  lines.  If  the  fetch  is  quite  small,  say  containing  only  one  or  two 
grid  points,  and  has  sharp  sides,  this  assumption  will,  thus,  cause  forecasts 
for  some  grid  points  to  be  much  too  high  and  for  other  grid  points  much  too 
low.  The  total  energy  in  the  entire  area,  however,  should  be  correct.  The 
only  way  that  this  total  energy  can  be  decreased  is  through  some  type  of  dissi- 
pation. Because  such  small  sharp  edged  fetches  seldom  occur  and  do  not  seem 
to  have  occurred  in  either  of  the  cases  studied,  this  assumption  does  not  appear 
to  be  at  fault.  The  edges  of  the  fetch  are  poorly  defined  so  that  all  the  sur- 
rounding grid  points  should  be  trading  part  of  their  energy  with  each  other  if 
the  model  allowed  it.  These  trades  should  not  change  the  total  energy  at  any 
one  grid  point  by  a large  percentage.  If  this  were  the  problem,  the  final 
significant  height  field  would  be  quite  irregular;  instead,  the  field  is  very 
easy  to  contour  with  smooth  maxima. 
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A second  dispersion  error  is  caused  by  not  having  forecast  low  enough 
frequencies  to  be  present.  These  lower  frequencies  travel  more  rapidly, 
thus  passing  the  observing  station  faster.  Along  with  this  is  the  fact  that 
the  model  propagates  waves  in  accordance  with  linear  wave  theory,  which  may 
be  slightly  too  low.  As  a more  complete  test  of  the  effects  of  dispersion, 
two  test  forecasts  were  prepared  for,.G6se  II.  In  the  first  of  these,  the 
assumption  was  made  that  the  wind  was  calm  after  step  18  (00Z,  December  18). 

A comparison  of  these  results"  with  the  previous  forecast  is  shown  in  Table  5-1 
for  grid  point  73.  This  comparison  shows  that  the  changing  wind  direction  was 
responsible  for  only  a small  part  of  the  error.  In  the  second  test,  besides 
the  same  assumption  of  a calm  wind  after  step  18,  the  frequencies  in  step  17 
were  decreased  by  l/4  of  their  value.  This  is  a slight  underestimate  of  the 
frequency  shown  in  Figure  4 -8(c)  so  that  it  results  in  a slight  overestimate 
of  the  propagation  velocity.  These  results  are  also  shown  in  Table  5-1  and 
again  account  for  another  small  part  of  the  error.  In  both  of  these  test 
forecasts  the  significant  heights  near  grid  point  73  decreased  fairly  evenly 
and  regularly.  Effectively,  therefore,  when  high  waves  are  experienced  over 
a large  area  surrounding  the  location  of  interest,  waves  cannot  propagate  out 
of  the  region  fast  enough  to  cause  the  extremely  rapid  decays  observed,  al- 
though reasonably  rapid  decays,  such  as  shown  for  ship  J in  Figure  4-5,  can 
occur.  One  further  test  should  be  made  in  a final  attempt  to  explain  these 
errors  in  terms  of  dispersion.  This  would  be  to  artificially  decrease  the  fre- 
quencies throughout  the  entire  computation.  Until  this  is  done  dissipation 
cannot  be  proved  conclusively. 

Since  the  observed  decay  may  not  completely  be  explained  by  dispersion  errors, 
the  need  for  a dissipation  function  must  be  considered.  As  stated  previously, 
the  Neumann  spectrum  was  not  devised  for  this  particular  application.  It  already 
considers  this  dissipation  term  for  the  area  of  generation. 

A simple  dissipation  function  could  be  added  very  easily  by  increasing  the  tab- 
ulated value 8 of  the  growth  function  and  decreasing  the  various  components.  It 
should  be  noted  that  some  dissipation  function  could  be  found  which  would  re- 
duce the  errors  in  Case  II  to  negligibility. 
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Table  5-1  Forecasts  of  Significant  Height  in  feet  for 
Grid  Point  73  Assuming  Special  Conditions 


Time 

Step 

Observed 

Regular  Forecast 

Calm  Wind ( 1 ) 

Calm  Wind^1)  (2) 
and  Fast  Propagation 

18 

35-3 

43.6 

43.6 

42.1 

19 

28.8(3) 

42.5 

41.5 

40.8 

20 

23. 9^ 

41.8 

40.7 

40.7 

21 

20.6 

43.0 

42.0 

39-9 

22 

19-5(3) 

42.6 

40.5 

38.7 

23 

16.1(3) 

41.7 

39-3 

38.7 

24 

13.6 

41.5 

38.9 

37.8 

(^All  winds  after  time  step  18  are  treated  as  calm. 

(2) a11  frequencies  decreased  by  25$  after  time  step  17  to  match  forecast 
errors  and  increase  propagation  velocity. 

(3) 

' Interpolated 


5-3  OUTLOOK 

One  of  the  major  objectives  of  the  present  paper  is  to  define  a numerical 
prediction  model  that  will  aid  the  future  development  of  wind-wave  forecasting 
through  objectivity.  It  should  be  relatively  easy  to  try  other  possible  spectra 
besides  the  one  used  and  conduct  an  extensive  verification  by  comparison  with 
observations.  Addition  of  the  effect  of  the  sea-air  temperature  difference  and 
other  stability  and  atmospheric  turbulence  criteria  on  generating  waves  should 
be  a worthwhile  improvement.  Dissipation  functions  can  be  easily  tested. 
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The  problem  of  specifying  the  initial  conditions  is  quite  serious  for  short 
range  forecasts.  For  example  as  shown  in  Figure  4-7,  it  took  a day  of  fore-  j 
casting  to  catch  up  with  an  initial  12  foot  error  in  significant  height.  The  / 
energy  that  could  not  have  come  from  a single  particular  wind  speed  should  not 
be  dropped.  It  might  help  if  a check  were  made  on  the  winds  at  surrounding 
grid  points  and  earlier  times  and  using  one  of  these  values  to  improve  the 
initial  estimate.  Another  simpler  and  possibly  more  promising  solution  might 

} 

be  to  estimate  the  principal  direction  of  the  initial  sea  and  use  this  with  / 
the  lowest  wind  speed  which  would  generate  a fully  developed  spectrum  having 
the  estimated  significant  height.  If  this  were  tried,  no  initial  wind  condi- 
tions would  be  needed.  / 

The  various  parameters  of  the  model,  such  as  grid  size  and  shape,  time  step 
length,  and  the  method  of  specifying  the  spectrum,  should  be  studied  in  some 
detail.  Perhaps  the  grid  size  should  be  increased  (more  space  between  points) 
and  the  extra  memory  space  thereby  saved  used  to  provide  two  tags  for  each 
component  so  that  the  exact  distance  past  and  to  the  side  of  each  grid  point 
could  be  specified.  This  system  might  allow  more  accuracy  in  location.  How- 
ever, this  increase  in  accuracy  will  be  limited  unless  the  curvature  of  the 
earth  is  also  considered.  As  suggested  previously,  other  methods  of  speci- 
fying the  spectra  should  be  considered. 

Most  of  the  possible  practical  applications  of  the  model  include  either  hind- 
cast  statistics  or  forecasting.  For  these  purposes,  a self  correcting  system 
should  be  added  so  that  all  wave  observations  can  be  added  as  available.  Other 
more  complete  methods  of  specifying  the  final  forecast  are  needed.  However, 
ship  motions  and  other  similar  parameters  which  depend  on  the  spectrum  can  be 
programmed  so  that  the  final  output  is  the  desired  parameter  directly. 
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I DECK  STACKING 

Monitor  Cards 
Program 
Data  Cards 


1 


Title  Card: 


APPENDIX  B PROGRAM 


Case  number  ending  in  Col.  8 


2 

3 

k 

5 . thru  64 
65  thru  X 

X thru  X + 28 

X + 29 

X + 30  thru  XX 


Print  Interval:  Interval  for  print-out  time  in  Cols.  2 - 12*. 
End  Card:  Number  of  final  step  in  Col.  2 - 12*,  number  of 
special  cards  ending  in  Col.  15 . 

Time : First  time  step  number 

Wind  Field:  Card  format  specified  below 

Special  Cards  (if  specified  in  Card  3)  time  step  in  Col  1-8*, 
Grid  point  number  in  Col.  14  - 16 

Initial  wave  fields  (if  time  = 0) : Card  format  specified  below 
Time 

Wind  Field  repeat  for  a11  time  steps 


II  CARD  FORMATS 


OOOOOOOQOOGOQQQQQQCQQCQOOQOOOOOOOOOOOOOOOOCOOOOOOOOOOOOOODOOOOOQOOOOOOGOOOOOOOOO 

I : 3 4 5 • 7 • 0 1C  IM2  13  14  H 16  »7  10  '9  23  ?!  22  23  24  ?S  20  27  20  24  30  31  32  33  34  3*  J6  37  30  39  40  4!  42  4 1 44  45  46  4?  40  4»y»  SI  *7  53  54  55  56  57  50  59  60  6!  62  61  M 65  66  67  60  69  70  71  72  ?3  ;t  75  7$  ?7  71  74  tO 

V V • D O D • V V • ODD*  VV*OOD.  VV*DoO«  VV*  ODD*  VV*ODD«  VV  • DOO»  VV»OOD«  VV*DOO<* 

»4  Ntl  N+2  N + 3 N+4  N+8  N+«  N4-7  N+8 

CARD  A.  WIND  FIELDS 


ooooooooooooooooooooooooooooooncoooooooooooooQoooooooooooooooooooooaoocaoGcooooo 

1 2 3 4 3 • ? • t 10  It  12  13  14  IS  It  77  It  ?3  7i  22  73  24  25  26  2?  20  23  3C  31  12  33  34  35  36  37  30  39  « 41  42  43  44  45  46  4\4I  49  M 51  52  53  54  55  50  57  SI  51  60  61  62  63  64  65  60  67  Cl  Cl  70  71  72  73  74  75  7*  77  71  79  00 
HH  • HH  • HH  • HH*  HH  • HH  • HH  • HH  « HH  . MM  • HH  . HH  HH  • HH  • HH  • HH  • HH  • HH  • 

N N+1  N+2  N+l  N+4  N+a  N+4  N+7  N+e  N + tO  N4-3 1 V0  + 12  N4-1S  H-M4  N4-5B  N-M7 

# 

\ 

CARD  B.  INITIAL  WAVE  FIELD 

where : 

w.  = wind  speed  in  \ 

odd.  = wind  direction  in  degrees  l at  location  n+m 

hh • = significant  wave  height  j 

n+m4  519  = relative  grid  point  number  over  water. 

*Use  decimal  point. 
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Table  B-l  Program 


FORTRAN 

2BAE0010 

MAIN  PROGRAM.  L.  BAER. WAVE  FCSTG2 

2BAE0020 

COMMON  SPACE. B I6FLD, E , STORGE . tMAX.O I R , FREQ. DELTFCJ, 

2BAE0030 

2TABCOL. TABROW. TABLE. JUMP, TIME. M I NDIS.NEVNTB, 

2BAE0040 

2NOOTAB 

2BAE0050 

DIMENSION  SPACE ( 86 A ) .BIGFLDI 86 A I.EI10.12 ) ,ENAX < 10 . 12 ) , 

2 BAE0060 

2STORGE< 20760) .01 R ( 12  I .FREQ ( 10 ) .DELTFQ ( 10 ) . TA0COL ( 25 1 . 

2BAE0070 

2TABROW1 1 1 ) .TABLE  12 5. 11)  .WINDDRI864)  .W1NDSPIB64)  • 

20AEOO8O 

2MINDISI 12 ) .NEVNTB 1 12 ) .NODTA01 12 1 

2BAE0090 

DIMENSION  TIMELI 1001 .LOCGRDIIOO) .FIRST (S19) 

28AE0100 

EOU I VALENCE (WINDDR. SPACE) , IWINDSP.BIGFLD ) 

2BAE0H0 

DIR< 1 ) «0, 

2BAE0120 

DIRI2 )*30. 

2BAE0130 

DIR ( 3 )*60. 

2BAE0140 

DIR ( A ) *90. 

2BAE0150 

DIR ( 5 ) *1 20. 

2BAE0160 

DIR ( 6 ) *150. 

2BAE0170 

DIR(7)»180. 

2BAE0180 

0IR(8 )*210. 

2BAE0190 

OIR 1 9 ) *240. 

2BAE0200 

DIR  1 10 ) *270. 

2BAE0210 

DIRI 11 1*300. 

2SAE0220 

DIR ( 12 ) *330. 

2BAE0230 

FREO ( 1 ) *.35 

2BAE0240 

FREOI 2)*1./5.31 

2BAE0250 

FREOI 3 ) *1. 27.77 

2BAE0260 

FREQI4 )«1 ./ 10.2  3 

28AE0270 

FRE0I5)*1./12.69 

2BAE0280 

FREOI 6 )*1./15.15 

2BAE0290 

FREOI7)*l./17.61 

2 8 AE 03 00 

FREO (8)«l./20.07 

2BAE0310 

FREO(9)*l./22.53 

2BAE0320 

FREOI 10)*. 04 

2BAE0330 

DELTFQ! 1 ) *.20 

2BAE0340 

DELTFO (2)*.5*(FREQ(1) -FREQ (311 

2BAE0350 

DELTFQ ( 3 ) * .5* I FREO 1 2 ) -FREQ (4 ) ) 

2BAE0360 

DELTFQ I 4)*«5*(FREQ(3)-FREQ(5) ) 

2BAE0370 

DELTFQ (3)*.3*(FREOI4)-FREQ(6) ) 

2BAE0380 

DELTFQ(6)*.5*(FREQ(5) -FREQ (7)1 

2BAE0390 

DELTFO(7)*.5*(FREQI6)-FREQ(8) ) 

2BAE0400 

DELTFQ(8)*.5*(FREQ(7)-FREQ(9)  ) 

2BAE04 10 

DELTFQI9 ) * .5* (FREOI 8 l-FREQ (10) ) 

2BAE0420 

DELTFQ ( 1 0 ) * .5* ( FREQ ( 9 ) - ( 1 . /27.46  ) ) 

2BAE0430 

TABCOLI1 )»0. 

2BAE0440 

TABCOL ( 2 ) *1 • 

2BAE0450 

T ABCOL (31*5. 

2BAE0460 

TABCOL ( 4 ) * 10. 

2BAE0470 

TABCOL (51*15. 

2BAE0480 

T ABCOL 1 6 ) *20. 

2BAE0490 

TABCOL (7) *25. 

2BAE0500 

T ABCOL (81*30. 

2BAE0510 

TABCOL (91*35. 

2BAE0520 

TABCOL ( 10 ) *40. 

2BAE0530 

i 

1 
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2BAE0540 

2BAE0550 

2BAE0560 

2BAE0570 

2BAE0580 

2BAE0590 

2BAE0600 

2BAE0610 

2BAE0620 

2BAEQ630 

2BAE0640 

2BAE0650 

2BAE0660 

2BAE0670 

2BAE0680 

2BAE0690 

2BAE0700 

2BAE0710 

2BAE0720 

2BAE0730 

2BAE0740 

2BAE0750 

2BAE0760 

2BAE0770 

2BAE0780 

2BAE0790 

28AE0800 

2BAE0810 

2BAE0820 

2BAE0830 

2BAE0840 

2BAE0850 

2BAE0860 

2BAE0870 

2BAE0880 

2BAE0B90 

2BAE0900 

2BAE0910 

2BAE0920 

2BAE0930 

2BAE094G 

2BAE0950 

2BAE0960 

2BAE0970 

2BAE0980 

2BAE0990 

2BAE1000 

2BAE1010 

2BAE1020 

2BAE1030 

2BAE1040 

2BAE1050 

2BAE1060 

2BAE1070 


TABCOLt 11 1*45. 

T ABCOL (121=50. 
TABCOL (131*60. 
TABCOLI 14) *70. 
TABCOL <J5) *80. 
TABCOL 1 16) *90. 
TABCOL (171*100. 
TABCOL (181*115. 
T ABCOL ( 1 9 ) * 130 . 
TABCOL (20 1*145. 
TABCOL (21 1*160. 
TABCOL (221*175. 
TABCOL (231=200. 
TABCOL (241*250. 
TABCOL ( 25 ) =400. 
TABROWI 1 1=16.88 
TABROVM  2 1=  23.6 
TABROWI 31=30.4 
T A BROW (4  1=  37.1 
T ARROW (51=43.9 
T ABROW (61=50.6 
TABROWI 7 1 =57.4 
TABROW ( 8 1 =64. 1 
TABROW (91=74.3 
TABROWI 101*84.4 
TABROWI 11 1=94.5 
DO  1 N=1 .25 
DO  1 L = 1 *11 
1 TABLE (N»L 1 =0. 
TABLE!  1.1 >*.l 
T ABLE  ( 1 . 2 1 * .1  5 
T ABLE ( 1 » 3 1 * . 2 
TABLF ( 1*4) =.35 
T ABLE ( 1 • 5 1 * .5 
T ABLE ( 1 .6  1 = .7 
T ABLE (1.71  = 1. 
TABLE (1*81*1.5 
TABLE ( 1.91*2.0 
TABLEIl. 101*3.1 
TABLF(1.11>=4.5 
TABLEI2.2  I-.6 
T ABL  E(2.3)*1.7 
T ABLE (2.41  = 1.8 
TABLEI2. 51*2.0 
TABLEI2. 61=2.5 
TABLE(2.7)=3. 

T ABLE (2.81*3.5 
TABLE(2. 91=4.0 
TABLE ( 2 « 1 0 1 =5 • 3 
TABLEI2.il  1=7. 

T ABLE (3.41=3.2 
TABLE ( 3 * 5 1 = 3 .4 
TABLF(3.6)*3.8 
TABLE(3,7)=4. 
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Program  (Continued) 


TABLE 

(3,81*4.5 

2BAE1080 

TABLE 

(3,9)«6. 

2BAE1090 

TABLE 

(3,101*8.7 

2BAEU00 

TABLE 

(3,111*12. 

2BAE1110 

TABLE 

(4,41*2. 

2BAE1120 

TABLE 

(4.51*4.4 

2BAE1130 

TABLE 

(4,61*4.5 

2BAE1140 

TABLE 

(4,71*4.5 

2BAE1150 

table 

(4,81*5. 

2BAE1160 

table 

(4,91*7.5 

2BAE1170 

table 

(4,101*11. 

2BAE1180 

TABLE 

(4,11 1=15. 

2BAE1190 

TABLE 

(5,51=5. 

2BAE1200 

TABLE 

15,61*5. 

2BAE1210 

TABLE 

(5.71*5. 

2BAE1220 

TABLE 

(5,81=5.4 

2BAE1230 

TABLE 

(5,91=8. 

2BAE1240 

table 

(5,101=12. 

2BAE1250 

TABLE 

(5,111*17. 

2BAE1260 

TABLE 

<6,51=5. 

2BAE1270 

TABLE 

(6,61*5.5 

2BAE1280 

TABLE 

(6,71=6. 

2BAE1290 

TABLE 

(6,81=6.4 

2BAE1300 

TABLE 

(6,91=9. 

2BAE1310 

TABLE 

(6.101*13. 

2BAE1320 

TABLE 

(6.11 1=20. 

2BAE1330 

TABLE 

(7,51=4.5 

2BAE1340 

TABLE 

(7,61*5.9 

2BAE1350 

TABLE 

(7,71=6.4 

2BAE1360 

TABLE 

(7,8 1 =6.9 

2BAE1370 

TABLE 

(7,91=9.4 

2BAE1380 

TABLE 

(7,101=14. 

2BAE1390 

TABLE 

(7,111*23. 

2BAE1400 

TABLE 

(8,61=6.6 

2BAE1410 

TABLE 

(8,71*7. 

2BAE1420 

TABLE 

[8,8 1=8. 

2BAE1430 

TABLE 

(8,91=10. 

2BAE1440 

TABLE 

(8,101*15. 

2BAE1450 

TABLE 

(8,111*25. 

2BAE1460 

tablf 

(9,61=7.7 

2BAE1470 

table 

(9,71=8. 

2BAE1480 

table 

(9,81*9. 

2BAE1490 

TABLE 

(9,91=10.4 

2BAE1500 

TABLE 

(9,101=16. 

2BAE1510 

table 

(9,111*25. 

2BAE1520 

tablf 

10,61=8.7 

2BAE1530 

TABLE 

( 10.71*9.1 

2BAE1540 

TABLE 

(10,81=9.8 

2BAE1550 

TABLE 

( 10,91*10.7 

2BAE1560 

TABLE 

( 10,101*17. 

2BAE1570 

TABLE 

(10,111=25. 

2BAE1580 

TABLE 

(11,61=9. 

2BAE1590 

TABLE 

(11,71=10. 

2BAE1600 

TABLE 

(11.81 *10.5 

2BAE1610 

4 ► 
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Table  B-l  Program  (Continued) 


TABLE (11,91=11. 

2BAE1620 

TABLE (11,101=17.5 

2BAE1630 

TABLEI1I. 111=25. 

2BAE1640 

TABLE! 12,61=7. 

2BAE1650 

TABLE (12,71=11. 

2BAE1660 

TABLEI12. 81=11. 

2BAE1670 

TABLE (12,91=11.4 

2BAE1680 

TABLEI12, 101=18. 

2BAE1690 

TABLEI12, 111=25. 

2BAE1700 

TABLEI13. 71=12. 

2BAE1710 

TABLE! 13.8 1=12. 

2BAE1720 

TABLE (13*91=12. 

2BAE1730 

TABLE! 13,101=18.5 

2BAE1740 

T ABLE (13,11 1=25. 

2BAE1750 

TABLEI14. 71=12. 

2BAE1760 

TABLF ( 14,8 1=1 2.4 

2BAE1770 

TABLE! 14,91=12.4 

2BAEI780 

TABLE! 14,101=19. 

2BAE1790 

TABLEI14. 111=25. 

2BAE1800 

TABLE! 15,71=12. 

2BAE1810 

TABLE(15,8)=12.6 

2BAE1820 

TA8LE(15»9)*12.8 

2BAE1830 

TABLEI15, 101=20. 

2BAE1840 

TABLEI15, 111  = 25. 

2BAE1850 

TABLEI16, 7)»11. 

2BAE1860 

TABLE(16,8)=12.9 

2BAE1870 

T ABLF (16,91=13.2 

2BAE1880 

TABLEH6*  101  = 21. 

2BAE1890 

T ABLE (16,11 1 = 25. 

2BAE1900 

TABLE (17,71=10. 

2BAE1910 

TABLE!  17,8  1=13.2 

28AE1920 

TABLE(17,9)=13.8 

2BAE 1930 

TABLEI17. 101=22. 

2BAE1940 

TABLE (17,11  1=25. 

2BAE1950 

TABLE(18,8)=13.5 

2BAE1960 

TABLE(18»9I=14.2 

2BAE1970 

TABLEI18, 101=23. 

2BAE1980 

TABLE!  18,11  1 = 25. 

2BAE1990 

TABLEI19, 81-14. 

2BAE2000 

TABLF! 19.9 1=15. 

2BAE2010 

TA8LE(19. 101  = 24. 

2BAE2020 

TABLE (19,111=25. 

2BAE2030 

TABLEI20, 8)=lt. 

2BAE2040 

TABLEI20, 91=15.8 

2BAE2050 

TABLE(20»101»25. 

2BAE2060 

TABLE(20.11  1 = 25. 

2BAE2070 

TABLE ( 21  ,81  = 13. 

2BAE2080 

TABLE (21  ,91=17. 

2BAE2090 

TABLEI21, 101=25. 

2BAE2100 

TABLE(21.11 1 = 25. 

2BAE2110 

TABLE (22,81=10. 

2BAE2120 

TABLE(22, 91=20. 

2BAE2130 

TABLE(22. 101=25. 

2BAE2140 

TABLEI22. 11  1*25. 

2BAE2150 
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Table  B-l  Program  (Continued) 


TABLE (23.91-23. 

2BAE2160 

T ABLE I 23 • 1 0 ) -25 

2BAE2170 

TABLEI23.il 1*25. 

2BAE2180 

TABLE  124.91-23. 

2BAE2190 

TABLE  (24.10L23. 

2BAE2200 

TABLfl24.il >«23. 

2BAE2210 

T ABL  F(23»9I»23. 

2BAE2220 

TABt  E(  23. 10). 25. 

2BAE2230 

TABLF(23,n  > -23 . 

2BAE2240 

MINOISIl  )>120 

2BAE2250 

MINDIS(2)-207.5 

2BAE2260 

MtNOISMI-207.5 

2BAE2270 

M 1 NO  1 $ 1 4 )• 1 20 

2BAE2280 

Ml  NO  13(31-207. 3 

2BAE2290 

Ml NO 1SI6I* 207. 3 

2BAE2300 

HINDIS! 71-120 

2BAE2310 

MI  NO  I SIB  1-207.5 

2BAE2320 

MIN0IS(9I«207.5 

2BAE2330 

HINDIS ( 1 0 ) -120 

2BAE2340 

MINDISdl  1-207.5 

2BAE2350 

MIN0IS(12)-207.5 

2BAE2360 

NEVNTB! 1 1-32 

2BAE2370 

NEVNTB ( 2 ) -31 

2BAE2360 

NEVNTBI3I  — 1 

2BAE2390 

NFVNTBI4I— 1 

2BAE2400 

NEVNTB ( 5 ) *-33 

2BAE2410 

NEVNTB ( 6 ) -“32 

2BAE2420 

NEVNTB  (7)— 32 

2BAE2430 

NEVNTB ( 8 1 =-3 1 

2BAE2440 

NE VNTB ( 9 ) * 1 

2BAE2450 

NEVNTB (101-1 

2BAE2460 

NEVNTB ( 1 1 1-33 

2BAE2470 

NFVNTB ( 12 ) -32 

2BAE2480 

NODTABI 1 )-32 

2BAE  2490 

NODTAB (21-32 

2BAE2500 

NODTABI 3 ) =3 1 

2BAE2510 

NODTAB ( 4 1 = — 1 

2BAE2520 

NODTABI 3 ) *-l 

2BAE2530 

NODTAB 1 6 1 *-33 

2BAE2540 

NODTAB  I 7 I— 32 

2BAE2550 

NODTAB I 8 ) =*32 

2BAE2560 

NODTAB I 9 ) -—31 

2BAE2570 

NODTABIlOI-l 

2BAE2580 

NODTABIll >-l 

2BAE2590 

NODTAB (121*33 

2BAE2600 

READ  INPUT  TAPE  5 . 1000 .N T I TLF . TMPRNT .ENDTME .NSPECT .NT I ME . 

2BAE2610 

l(WlNDSPIN) .WINDDRIN) .N-l .519) 

2BAE2620 

1000 

FORMAT  I I8/F12.0/F12.0. I3/I8/I 18F4.0)  ) 

2BAE2630 

IF  (NSPFCTI4000. 4000. 4001 

2BAE2640 

4001 

READ  INPUT  TAPE  5 , 4002 • I T I MEL  I K) .LOCGRD I K ) .K- 1 .NSPECT ) 

2BAE2650 

4002 

FORMAT  I F8 .0 . 1 8 I 

2BAE2660 

C 

TO  INSPECT  COMPLETE  SPECTRA.  PUT  NUMBER  OF 

2BAE2670 

C 

POINTS  TO  BE  COVERED  AS  NSPECT  AND  ADD  A CARD 

2BAE2680 

c 

AFTER  FIRST  WINDFIELD  FOR  EACH  POINT.  MAX 

2BAE2690 
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NUMBER  OF  POINTS=100.  PRINTOUT  TIME  ONLY 

2BAE2700 

4000 

CALL  SDH  131) 

2BAE2710 

REWIND  31 

2BAE2720 

IF  INTIME )2.3»2 

2BAE2730 

3 

READ  INPUT  TAPE  5 .2000 . ( F I RST ( N ) .N  = 1 . 5 19 ) 

2BAE2  740 

2000 

FORMAT! 18F4.1 ) 

2BAE2750 

TIME=0. 

2BAE2760 

DO  2020  1=1.20760 

2BAE2770 

2020 

STORGE  ( I ) = 0. 

2BAE2780 

DO  2009  NGR ID= 1 1 5 1 9 

2BAE2790 

IF  I WINDSP INGRID) -10. 12009. 20 10. 2010 

2BAE2800 

2010 

IF  (WINDSP I NGR ID) -56. 12040,2040.2041 

2BAE2810 

2 041 

WRITE  OUTPUT  TAPE  6 , 502 . T I ME .NGR I D.W I NDSP I NGR I D ) 

2 8AE2820 

W I NDSP I NGRID)=56. 

2BAE2830 

2040 

CALL  SPREAD  I NGR I D ,NGR I D ,0 ) 

2BAE2840 

WI NDSP (NGR ID) =1.688944*WINDSP INGRID) 

2BAE2850 

CALL  COMPEl WINDSP (NGR ID) , W INDDR (NGR I D ) • F IRST ( NGR I D 1 I 

2BAE2860 

CALL  REPACK  ( NGR I D .NGR ID ,0 ) 

2BAE2870 

2009 

CONTINUE 

2BAE2660 

GO  TO  2050 

2BAE2890 

2 

PRINT  500 

2BAE2900 

500 

FORMAT ( 33H1  PLEASE  MOUNT  RESERVE  TAPE  ON  Cl) 

2BAE2910 

PAUSE 

2BAE2920 

REWIND  31 

2BAE2930 

101 

DO  6 N=1 .NT IME 

2BAE2940 

6 

READ  TAPE  3 1 . T I ME .STORGE 

2BAE2950 

100 

TIME®  TIME  + 1. 

28AE2960 

IF (TIME-FLOATFINTI ME) 17,8.7 

2BAE2970 

7 

CALL  POUMP 

2BAE2980 

GO  TO  25 

2BAE2990 

S 

DO  9 NGR  I D = 1 ,519 

2BAE3000 

I F ( W I NDSP ( NGR I D ) -1 0. ) 9,10.10 

2BAE3010 

10 

I F ( W I NDSP ( NGR I D ) -56. ) 40.  40.41 

2BAE3020 

41 

WRITE  OUTPUT  TAPE  6 . 502 » T I ME . NGR ID ,W I NDSP ( NGR I D ) 

2BAE3030 

502 

FORMAT (1HO/8HOT1ME  = F 12 .0 . 10HGR ID  PT  = 1 8 • 8HSPEED  = F12.2) 

2BAE3040 

WINDSP! NGRID)=56. 

2BAE3050 

40 

CALL  SPREAD  I NGR I D ,NGRI D ,0 ) 

2BAE3060 

WINDSP(NGRID)=1.688944*WINDSP(NGRID) 

2BAE3070 

CALL  COMPE (WINDSP ( NGR 1 D ) ,W INDDR ( NGR I D ) . 0 ) 

2BAE3060 

CALL  REPACK  INGRID,  NGR ID,  0) 

2BAE3090 

9 

CONTINUE 

2BAE3100 

CALL  PRPAGT 

2BAE3110 

2050 

WRITE  TAPE  31. TIME, STORGE 

2BAE3120 

I F ( MODE ( T I ME  » TMPRNT ) ) 13, 12. 13 

2BAE3130 

PRINT  ROUTINE  BEGINS  AT  LOCATION  12 

2BAE3140 

12 

DO  14  N=l,864 

2BAE3150 

14 

SPACE  ( N ) =0. 

20AE316O 

DO  15  N=1 .12 

28AE3170 

DO  15  M= 1*10 

2BAE3180 

CALL  SPREADIM.N.l ) 

28AE3190 

DO  15  L= 1 ,864 

2BAE3200 

15 

SPACE (L ) =SPACE (L ) +BIGFLD(L ) 

2BAE3210 

DO  16  L*1 ,864 

2BAE3220 

16 

SPACE ( L ) = 2.83*SORTF (SPACE! L I ) 

28AE3230 

B-8 
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DO  17  M* 1 *17*16 
IFIM-l 119,19,18 

19  WRITE  OUTPUT  TAPE  6.20.NTITLE 

20  FORMAT ( 34H1BAER  WAVE  FORECAST  LEFT  SECT ION/6H  CASE  lb) 

GO  TO  21 

18  WRITE  OUTPUT  TAPE  6.22 

22  FORMAT ( 33H1 BAER  WAVE  FORCAST  RIGHT  SECTION) 

21  WRITE  OUTPUT  TAPE  6 ,2 3 . T IME , ! SP ACE < N ) , SPACE < N+l  I , 

1 SPACE ( N+2 ) .SPACE  I N+3) . SPACE ( N+4 1 , SPACE ( N+S ) .SPACE (N+6)  . 
1SPACEIN+7) .SPACE (N+8), SPACE (N+9) .SPACE (N+10) .SPACE (N+ll  I , 

1 SPACE (N+12 ) .SPACE ( N+13 ) .SPACE (N+ 14) , SPACE (N+ 15 ) ,N=M.864, 

132) 

23  FORMAT tIHOlOX.llHTIME  STEP  F 12 .0/ ( 1H0 /3X . 1 6F7 . 1 ) ) 

17  CONTINUE 

IF  (NSPECTI13.13 .4009 

4009  DO  4010  K=1,NSPECT 

IF  IT1MEL(K)-T1ME)4010. 4012. 4010 

4012  CALL  SPREAD  I LOCGRD ( A ) , LOCGRO ( X) ,0 ) 

NGRID=LOCGRD(K ) 

WRITE  OUTPUT  TAPE  6 .401 3 . T IMF .NGRID. I ( E ( N ,M ) ,N=1 , 

110) ,M=1 .12  ) 

4013  FORMAT (8H1SPECIAL/8H0TIME  = F 1 5 . 0 . 1 2HGR I DPOI NT  = I8/I5E23.8)) 

4010  CONTINUE 

13  WRITE  OUTPUT  TAPE  6.24. TIME 

24  FORMAT ( 1HO/1HO/1HO/11HOTIME  STEP  F 12 .0 . 8HF I N I SHED ) 

IF  ( T I ME-ENOTME ) 30.25.25 

30  READ  INPUT  TAPE  5 . 3 1 .N T I ME . ( W I NDSP ( N I , W I NDDR ( N ) ,N 
1=1.510) 

31  FORMAT! 18/ ( 18F4.0 ) I 
GO  TO  100 

25  CALL  UNLOAD ( 3 1 I 
PRINT  501 

501  FORMAT  (32H1  PLEASE  RESERVE  TAPE  ON  UNIT  Cl) 

CALL  EXIT 
END 

* FORTRAN 

SUBROUTINE  COMPE ( W I NDSP , W I NDDR . F I RST ) 

C COMPE.  SUBROUTINE  TO  COMPUTE  ALL  120  NEW  E VALUES 

COMMON  SPACE. BIGFLD.E.STORGE.E MAX, DIR. FREQ. DELTFO, 

2TABC0L. TABROW, TABLE. JUMP.TIME. MIND  IS, NEVNT8, 

2N0DTAB 

DIMENSION  SPACE ( 864 ) .BIGFLD! 864 ),E< 10.121 .EMAX! 10.12) . 
2ST0RGF! 20760)  .D I R ( 1 2 ) , FRFO ( 10 ) .DEL TFQ ( 10), T ABCOL ( 25 ) , 

2 TABROW! 11  I .TABLE! 2 5.11  I ,DWN0DR(864) , DWNDSP l 864 ) , 

2MINDISI12) ,NE VNTB ( 1 2 ) .NODTAB112) 

EOU I VALENCE (DWNDDR. SPACE ) .(DWNDSP, BIGFLD) 

S I GMA  = 0 • 

DO  200  N* 1 , 10 
DO  200  M*1  . 12 
200  EMAX(N«M)=0. 

DO  210  JDIR=1.12 

SIGDR=ABSF(WINDDR-DIR( JDIR) 1-90. 

I F ( S I GDR 1202.202,205 

205  SIGDR=ARSF(WINDDR-360.-DIR( JDIR) )-90. 
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2BAE3240 
2BAE3250 
2BAE3260 
2BAE3270 
2BAE3260 
2BAE3290 
2BAE3300 
2BAE3310 
2BAE3320 
2BAE3330 
2BAE3340 
2BAE3350 
2BAE3360 
2BAE3370 
2BAE3380 
2BAE3390 
2BAE3400 
2BAE3410 
2BAE3420 
2BAE3430 
2BAE3440 
2BAE3450 
2BAE  3460 
2BAE3470 
2BAE3480 
2BAE3490 
2BAE3500 
2BAE3510 
2BAE3520 
2BAE3530 
2BAE3540 
2BAE3550 
2BAE3560 
2BAE3570 
2BAE3580 
2BAE3590 
2BAE3600 
2BAE3610 
2BAE3620 
2BAE3630 
2BAE3640 
2BAE3650 
2BAE3660 
2BAE3670 
2BAE3680 
2BAE3690 
2BAE3700 
28AE3710 
2BAE3720 
2BAE3730 
2BAE3740 
2BAE3750 
2BAE  3 760 
2BAE3770 


n. 
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! 


I 


f 


IFISIGDR1202.202.206 

206  SIGDR=AB5F(WINDDR+360.-DIR(JDlRI  1-90. 
IFISIGDRI202.202.210 

202  DO  204  1=1,10 

20 A S I GMA=S I GMA+F ( I « JD 1 R ) 

BETA=SIGDR+90. 

IF (BETA-75. 1215.215.216 
216  DELD I R=0. 523598 78 *( BETA-75. )/ 30. 

BETA=75.+(BFTA-75. )/2. 

GO  TO  203 

215  DELDIP=0. 52359878 

203  DO  201  JFREQ= 1,10 

211  EMAXI JFREQ.  JD1R)  = I 1 03 . *E XPF ( -2 .* ( 32 . 1 6 / ( FREQ ( JFREQ ) » 
16.283 185 3+WINDSP ) ) **2 ) *DEL TFQ < JFREQ ) *DELD I R/ ( 6 .283 1853 
1*FREQ( JFREQ) ) **6 ) * ( 1 . + ( . 5+0. B2+EXPF ( -0. 5» ( 

16.28  3185  3*FREQ( JFREQ ) *W I NDSP/ 32 . 16  I **4 ) I *COSF ( 0 .03490 7 
1*  BETA)+0.32*EXPF(-0.5*(6.2831853*FREQ( JFREQ) 

1 *W  I NDSP/32 . 16 ) *+4 ) *COSF ( 0.06981 3 + BET  A ) ) 

201  CONTINUE 
210  CONT1NUF 

IF  (FIRST  1300,300,301 

301  TABVAL* (FIRST**2 ) /8. 

GO  TO  302 

300  DO  221  KCOL=l ,25 

IF  I TABCOLIKCOL  > -SIGMA  1221, 222.222 

221  CONTINUF 

222  DO  223  KROW=l.ll 

I F ( TA3ROW (KROW)-WINDSP)223 .2 24.224 

223  CONTINUE 

224  COL1* (TABLE (KCOL.KROW-l ) -TABLE (KCOL-1. KROW-1 
1) )*(SIGMA-TABCOL(XCOL-l ) l/(TABCOL(KCOL)-TABCOL 
1 (KCOL-1 ) l+TARLE (KCOL-1 .KROW-1 ) 

COL2=( TABLE (KCOL.KROW) -TABLE ( XCOL-1 ,KROW ) )» 

1 (SIGMA-TABCOL (KCOL-1 ) ) / ( T ABCOL ( KCOL ) -T ABCOL 
1 (KCOL-1 ) )+TARLE(KCOL-l.KROW> 
TA8VAL=(COL2-COLl)*(WINDSP-TABROW( KROW-1 ) ) / 

2 ( T ABROW ( KROW ) -TABROW (KROW-1 ) ) +COL1 

302  DO  230  Ls 1 ,10 
DO  230  K=1 .12 

EMAXI L.K I ®EMAX(L,K)-E(L,K ) 

IF  (EMAX(L.KI)  231.230.230 
231  EMAX ( L »K I =0. 

230  CONTINUE 
SUM  = 0. 

C CHOOSE  DIR  NEAR  WINDDR 

DO  240  JDIRM.12 
IF (DIR(JDIR)-WINDDR 1 240.24 1,241 

240  CONTINUF 

241  IFIJDIR— 3)260.269,261 
261  IF(JDIR-10)262,268,263 
260  IF { JDIR-21264.265,265 
263  IFIJCIR-ll  1266.266.267 
266  KAMO 

KB*  12 


2BAE3780 
2BAE3790 
2BAE3800 
2BAE3810 
2 BAE3820 
2BAE  3830 
2BAE3840 
2BAE3850 
2BAE3860 
2BAE3870 
2BAE3880 
2BAE3890 
2BAE3900 
2BAE3910 
2BAE3920 
2BAE3930 
2BAE3940 
2BAE3950 
2BAE3960 
2BAE3970 
2BAE3980 
2BAE3990 
2BAE4000 
2BAE4010 
2BAE4020 
2BAE4030 
2BAE4040 
2BAE4050 
2BAE4060 
2BAE4070 
2BAE4080 
2BAE4090 
2BAE4100 
2BAE4110 
2BAE4120 
2BAE4130 
2BAE4140 
2BAE4150 
2BAE4160 
2BAE4 1 70 
2BAE4180 
2BAE4190 
2BAE4200 
2BAE4210 
2BAE4220 
2BAE4230 
2BAE4240 
2BAE4250 
2BAE4260 
2BAE4270 
2BAE4280 
2BAE4290 
2BAE4300 
2BAE43 10 


I 

i 
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KC*9 

2BAE4320 

ICO®  1 

2BAE4330 

KE  = 8 

2BAE4340 

KF®2 

2BAE4350 

GO  TO  270 

2BAE4360 

267 

KA  = 11 

2BAE4370 

KB*  1 

2BAE4380 

KC*  1 0 

2BAE4390 

K0*2 

2BAE4400 

KE  *9 

2BAE4410 

KF  = 3 

2BAE4420 

GO  TO  270 

2BAE4430 

262 

KA® JD I R— 1 

2BAE4440 

KB® JO  I R+ 1 

2BAE4450 

KC=JDIR-2 

2BAE4460 

K0= JD I R+2 

2BAE4470 

KE=JDIR-3 

2BAE4480 

KF=JDIR+3 

2BAE4490 

GO  TO  270 

2BAE4500 

266 

KA®  1 

2BAE4510 

KB  = 3 

2BAE4520 

KC*12 

2BAE4530 

KD=4 

2BAE4540 

KE*  1 1 

2BAE4550 

KF*5 

2BAE4560 

GO  TO  270 

2BAE4570 

264 

KA  = 12 

2BAE4580 

KB*2 

2BAE4590 

KC*11 

2BAE4600 

KD*3 

2BAE4610 

KE*10 

2BAE4620 

KF  *4 

2BAE4630 

GO  TO  270 

2BAE4640 

268 

KA*9 

2BAE4650 

KB®  1 1 

2BAE4660 

KC*8 

2BAE4670 

K0*12 

2BAE4680 

KE*7 

2BAE4600 

KF*1 

2BAE4700 

GO  TO  270 

2BAE4710 

269 

KA«2 

2BAE4720 

KB*4 

2BAE4730 

KC®1 

2BAE4740 

KD*5 

2BAE4750 

KE*12 

2BAE4760 

KF  *6 

2BAE4770 

270 

00  242  JFREQ=ltlO 

2BAE4780 

SUM=SUK+EMAX( JFREQ.JDIR) 

2BAE4790 

El JFREO.JDIR) =t 1 JFREQ.JDIR 1+EMAxi JFREQ. JDIR) 

2BAE4800 

IF  I SUM-TABVAU 243.245. 244 

2BAE4810 

243 

SUM=SUM+EMAX<  JFREQ.KA) 

2BAE4820 

E(JFRE0,KA)»E( JFREO  »K A ) + EMAX (JFREO.KA) 

2BAE4830 

IF (SUM-TABVAL 1246.245. 247 

2BAE4840 

246 

SUM*SUM+EMAX ( JFREO .KB  1 

2BAE4850 

C 
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E( JFREQ,KB)=E( JFREQ , KB ) +EMAX ( JFREQ, KB) 

IE ( SUM-TABVAL ) 248  *245*249 

248  SUM=SUM+EMAX (JFREQtKC) 

E (JFREQtKC )=E ( JFREQtKC )+EMAX( JFREQtKC ) 

IF ( SUM-TABVAL ) 250 t 245 t 251 

250  SUM=SUM+EMAX ( JFREQtKD) 

l E( JFREQ, KD)=E( JFREQ, KDl+EMAXI JFREQtKD) 

* I F ( SUM-T  ABVAL ) 252  t245t253 

252  SUM=SUM+FMAX (JFREQtKE) 

E( JFREQtKE)=F( JFREOtKE l+EMAXI JFREQtKE) 

IF  I SUM-TABVAL I 254 .245 1 25  5 

254  SUM-SUM+FMAX ( JFREQtKF) 

E I JFREQtKF I =E ( JFREQ.KF ) +EMAX I JFREQtKF I 
IF (SUM-T ABVAL I 242. 245, 25 7 

244  E I JFREQ. JDIRI=E( JFREQ.JDIR 1-1 SUM-TABVAL I 
GO  TO  245 

247  E (JFREQ, KA)=E I JFREQ.KAI- (SUM-TABVAL) 

GO  TO  245 

249  E( JFREQ. KB) =E( JFREQ, KB)- (SUM-TABVAL I 
GO  TO  245 

251  E(JFRFQ,KC)=F( JFREQtKC I- (SUM-T ABVAL I 
GO  TO  245 

253  E( JFREQ, KD)=E( JFREQ, KDI-( SUM-T ABVAL I 
GO  TO  245 

255  E (JFREQtKE) =E( JFREQtKE I- (SUM-TABVAL ) 

GO  TO  245 

257  E( JFREQtKF )=E( JFREQ.KF (-(SUM-TABVAL I 
GO  TO  245 
242  CONTINUE 

245  RETURN 
END 

* FORTRAN 

SUBROUTINF  PRPAGT 

COMMON  SP ACE, B IGFLD ,E tSTORGE tEMAX , DIR, FREQ. DELTFQ, 
2TABCOL.T  ABROW .TABLE. JUMP, TIMF .MIND  IS. NEVNTB, 

2NODTAB 

DIMENSION  SPACE (864) , 3 I GFLD ( 864 ) , E ( 10 . 12 ) tEMAX ( 10.12) . 
2STORGE (20760) .DIR ( 12) .FREQ ( 10) . DEL  TFQ ( 10 ) . TABCOL ( 2 5 ) t 
2TABROW1 11 ) .TABLE (25 .11) ,WINDDR(864I .WIN0SP(864) t 
2MINDISI 12 ) .NEVNTB ( 12  I .NODTABI 12 ) 

EQUI VALENCE (WINDDR. SPACE ) . (WINDSP.BIGFLD ) 

EVEN = MODE ( T IME .2 • ) 

DO  601  M0IR=1.12 
00  602  MFREQ= 1,10 

CALL  JUMP TM ( FREQ ( MFREO ) ,MI ND I S ( MD I R ) ) 

I F ( JUMP  1602,604,602 
604  CALL  SPREAD(MFREQ.MDIR.l) 

IF (EVEN  1609,608 .609 

608  LOCMOD= NEVNTB ( MD I R ) 

GO  TO  610 

609  LOCMOD=NODTAB(MDIR ) 

610  DO  605  MGR  I D= 1 ,864 
NUGRID=MGRID-LOCMOD 

SPACE ( MGR  I D l*BIGFLD( NUGR ID ) 


2BAE4860 

2BAE4870 

2BAE4880 

2BAE4890 

2BAE4900 

2BAE4910 

2BAE4920 

2BAE4930 

2BAE4940 

2BAE4950 

2BAE4960 

2BAE4970 

2BAE4980 

2BAE4990 

2BAE5000 

2BAE5010 

2BAE5020 

2BAE5030 

2BAE5040 

2BAE5050 

2BAE5060 

2BAE5070 

2BAE5080 

2BAE5090 

2BAE5100 

2BAE5110 

2BAE5120 

2BAE5130 

2BAE5140 

2BAE5150 

2BAE5160 

2BAE5170 

2BAE5180 

2BAE5190 

2BAE5200 

2BAE5210 

2BAE5220 

2BAE5230 

2BAE5240 

2BAE5250 

2BAE5260 

2BAE5270 

2BAE5280 

2BAE5290 

2BAE5300 

2BAE5310 

2BAE5320 

2BAE5330 

2BAE5340 

2BAE5350 

2BAE5360 

2BAE5370 

2BAE5380 

2BAE5390 
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605  CONTINUE 

CALL  REPACKIMFREQ.MDIR.l > 

602  CONTI NUF 
601  CONTINUF 
611  RETURN 
END 

FORTRAN 

SUBROUTINE  JUMPTM ( FREQCY  »MND I ST  I 
IT  IS  TIME  TO  JUMP  IF  JUMP«0 

COMMON  SPACE.BIGFLD.E.STORGE.EMAX.DIR.FREQ.DELTFQ, 
2TABCOL  . T ABROW. TABLE  .JUMP .TIME.MINDIS.NEVNTb. 

2NODTAB 

DIMENSION  SPACE (864) .BIGFLDI864) .E( 10.12) .EMAXI 10.12) . 
2STORGFI 20760) .D I R <12 ) .FREQ ( 10 ) .DELTFOI 10  I .TABCOL < 25  I . 
2TABROW (11) . TABLE <25.1 1 ) »W I NDDR (864) .W I NDSP ( 864 ) . 
2MINDIS<12) .NEVNTBI12) .NODTABI 12) 

EQU I VALENCE (WINDDR.SP ACE) . < W I NDSP .8 1 GFLD ) 

DISMIN'MNDIST 
TRAVD5=3 • /FREQCY 
DSLAST=(TI ME-1 • ) *TRAVDS 
PASTDS=MODF<DSLAST.DISMIN> 

IF (PASTDS+TRAVDS-DISMIN) 1001. 1002. 1002 


o 


28AE5400 
2BAE5410 
2BAE5420 
2BAE5430 
2BAE5440 
2BAE5450 
2BAE  5460 
2BAE5470 
2BAE54B0 
2BAE5490 
2BAE5500 
2BAE5510 
2BAE5520 
2BAE5530 
2BAE5540 
2BAE5550 
2BAE5560 
2BAE5570 
2BAE5580 
2BAE5590 
2BAE5600 
2BAE5610 


1001  JUMP= 1 

2BAE5620 

RETURN 

2BAE5630 

1002  JUMP=0 

2BAE5640 

RETURN 

2BAE5650 

END 

2BAE5660 

* FORTRAN 

2BAE5670 

FUNCTION  LANDSE ( MCNT ) 

2BAE5680 

MWHOLE* 

MCNT/36 

2BAE5690 

MPART  *XMODF (MCNT . 36 ) 

2BAE5700 

LANDSE= 

LSHIFTIMWHOLE. MPART  I 

2BAE5710 

RETURN 

2BAE5720 

END 

2BAE5730 

# FAP 

2BAE5740 

ENTRY 

SPREAD 

2BAE5750 

ENTRY 

REPACK 

2BAE5760 

ENTRY 

lshift 

2BAE5770 

SPREAD  SX A 

RETURN.  1 

SAVE  INDICES 

2BAE5780 

SX  A 

RETURN+1.2 

2BAE5790 

SX  A 

RETURN+2 .4 

2BAE5800 

CLA* 

1.4 

2BAE5810 

STO 

MFREQ 

MFREQ  IF  CODE  U MGR  ID  IF  CODEO 

2BAE5820 

CLA* 

2.4 

2BAE5830 

STO 

MDIR 

2BAE5840 

CLA* 

3.4 

TEST  CODE 

2BAE5850 

TN2 

C0DF1 

DO  BIG  FLD  ON  TRA 

2BAE5860 

CLA 

MFRFQ  MGR  ID 

2BAE5870 

SUB 

FID 

2BAE5880 

XCA 

2BAE5890 

MPY 

F40A 

2BAE5900 

XCA 

2BAE5910 

STD 

BASFO 

2BAE5920 

AXT 

0.1 

IR1  TO  40 

2BAE5930 

6 
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AXT 

0.2 

IR2  TO  120 

2BAE5940 

BASEO  TXI 

*+1.1 .0 

DECR  FROM  *-3 

2BAE5950 

SPR  CAL 

STORGE.l 

PACKED  WORD  TO  AC 

2BAE5960 

ANA 

MAStCL 

LEFT  WORD  IN  AC  LOGICAL 

2BAE5970 

ARS 

24 

TO  AC  ADDR 

2BAE5980 

TSX 

FLOAT. 4 

2BAE5990 

SUB 

OCTF 

2BAE6000 

STO 

E.2 

2BAE6010 

XCA 

2BAE6020 

FMP 

E . 2 

2BAE6030 

STO 

E.2 

2BAE6040 

CAL 

STORGE.l 

2BAE6050 

ANA 

MASKC 

CENTER  WORD 

2BAE6060 

ARS 

12 

2BAE6070 

TSX 

float. 4 

2BAE6080 

SUB 

OCTF 

2BAE6090 

STO 

E-1.2 

2BAE6100 

XCA 

2BAE6110 

FMP 

F-I.2 

2BAE6120 

STO 

F-1.2 

2BAE6130 

CAL 

STORGE.l 

2 BAE614C 

ANA 

MASKR 

RIGHT  WORD 

2BAE6150 

TSX 

FLOAT. 4 

2BAE6160 

SUB 

OCTF 

2BAE6170 

STO 

E-2,2 

2BAE6160 

XCA 

2BAE6190 

FMP 

E-2.2 

2BAE6200 

STO 

E-2.2 

2BAE6210 

TXI 

*+1,1.1 

2BAE6220 

TXI 

*+1,2.3 

2BAE6230 

TXL 

SPR .2,119 

2BAE6240 

RETURN  AXT 

0,1 

ADDR  FROM  SPREAD  OR  REPACK 

2BAE6250 

AXT 

0,2 

SPREAD+ 1 

2BAE626C 

AXT 

0.4 

SPREAD+2 

2BAE6270 

TRA 

4,4 

EXIT 

2BAE6280 

CODE  1 CLA 

MDIR 

SPREAD  864 

2BAE6290 

SUB 

FID 

2BAE6300 

XCA 

2BAE6310 

MPY 

F10A 

2BAE6320 

XCA 

2BAE6330 

ADD 

MFRFO 

2BAE6340 

SUB 

FID 

2BAE6350 

STZ 

WORK 

CHOOSE  MASK  ROUTINE 

2BAE6360 

STD 

WORK 

2BAE6370 

CLM 

2BAE6380 

LDO 

WORK 

2BAE6390 

LLS 

0 

2 BAE  6400 

DVP 

OCT3D 

REMAINDER  IN  AC  DECR 

2BAE6410 

STO 

WORK 

2BAE6420 

CLM 

2BAE6430 

LLS 

18 

2BAE6440 

XCA 

2BAE64S0 

STD 

BASE  1 

2BAE6460 

CLA 

WORK 

2BAE6470 

• 1 
rj 

I 
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ARS 

18 

2BAE6480 

TZE 

LMASK 

NO  REMAINDER,  CHOOSE  LEFT  WORD 

2BAE6490 

SUB 

0CT1A 

2BAE6500 

TNZ 

RMASK 

TRA  IF  REM=2,  CHOOSE  RIGHT  WORD 

2BAE6510 

CLA 

ADMSKC 

REM=1,  CHOOSE  CENTER  WORD 

2BAE6520 

STA 

MSKPL 

2BAE6530 

CLA 

AD12 

2BAE6540 

STA 

SHFTPL 

2BAE655C 

TRA 

THENSP 

* 

2BAE6560 

RM ASK  CLA 

ADMSKR 

2BAE6570 

STA 

MSKPL 

2BAE6580 

CLA 

ADO 

2BAE6390 

STA 

SKFTPL 

2BAE6600 

TRA 

THENSP 

2BAE6610 

LMASK  C. A 

ADMSKL 

2BAE6620 

STA 

MSKPL 

2BAE6630 

CLA 

AD24 

2BAE6640 

STA 

SHFTPL 

28 AE  6630 

THENSP  AXT 

0.1 

FOR  OUTPUT  TO  864 

2BAE6660 

AXT 

0.2 

FOR  INPUT  TO  20720 

2BAE6670 

BASE1  TXI 

*+1.2,0 

DECR  FROM  CODE1+13 

2BAE6680 

BIGSPR  SXD 

INDEX, 1 

2BAE6690 

TSX 

slandse.a 

2BAE6700 

TSX 

INDFX.O 

2BAE6710 

TNZ 

BIGSEA 

SEA  IF  1 IN  AC 

2BAE6720 

STZ 

BIGFLD.l 

LAND 

2BAE6730 

TRA 

BIGSEA  CAL 

TEST+1 

2BAE6740 

STORGE  » 2 

2BAE6750 

mskpl  ana 

0 

ADDR  OF  MASK  FROM  RMASK-2 , + 1 ,0R 

LMASK2BAE  >760 

SHFTPL  ARS 

0 

ADDR  FROM  RMASK+3  ETC 

2BAE6770 

TSX 

FLOAT .4 

2BAE6780 

SUB 

OCTF 

2BAE6790 

STO 

BIGFLD.l 

2BAE6800 

XC  A 

2BAE6810 

FMP 

BIGFLD.l 

2BAE6820 

STO 

BIGFLD.l 

2BAE6830 

| TEST  TXI 

*+1.2 . AO 

2BAE6840 

TXI 

*+1.1.1 

2BAE6850 

TXL 

BIGSPR. 1.863 

2BAE6860 

TRA 

RETURN 

2BAE6870 

REPACK  SX A 

RETURN.  1 

SAVE  INDICE! 

2BAE6880 

SXA 

RETuRN+1 .2 

2BAE6890 

SX  A 

RETuRN+2 .4 

2BAE6900 

CLA* 

1 .4 

2BAE6910 

STO 

MFREO 

MFREO  FOR  CODE  1 ,MGR I D FOR  CODEO 

2BAE6920 

<:la* 

2.4 

2BAE6930 

STO 

MDIR 

2BAE6940 

CLA* 

3.4 

TEST  CODE 

2BAE6950 

TNZ 

RCODE1 

DO  BIG  FLD  ON  TRA 

2BAE6960 

CLA 

MFRFQ 

MG*R  ID 

2BAE6970 

SUB 

FID 

2BAE6980 

XCA 

2BAE6990 

HP  Y 

F40A 

2BAE7000 

I 

2BAE7010 

£ 

■ 
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Table  B-l  Program  (Continued) 


STD 

RBASEO 

2BAE7020 

AXT 

0.1 

IR1 

TO 

40 

2BAE7030 

AXT 

0.2 

IR2 

TO 

120 

2BAE7040 

RBASEO 

TXI 

*+1.1.0 

DECR 

! FROM  *-3 

2BAE7050 

REP 

CLA 

E.2 

2BAE7060 

TSX 

SSQRT.4 

2BAE7070 

ADD 

OCTF 

2BAE7080 

TSX 

FIX. A 

BINARY  pt 

AFTER  BIT  5 

2BAE7090 

STO 

WORK 

2BAE7100 

CAL 

WORK 

2BAE71 10 

ANA 

MASKR 

2BAE7120 

ALS 

24 

2BAE7130 

SLW 

WORK 

2BAE7140 

CLA 

E-1.2 

2BAE7150 

TSX 

SSQRT.4 

2BAE7160 

ADD 

OCTF 

2BAE7170 

TSX 

FIX. 4 

2BAE7180 

STO 

WORK + I 

2BAE7190 

CAL 

WORK+1 

2BAE7200 

ANA 

MASKR 

2BAE7210 

ALS 

12 

2BAE7220 

SLW 

WORK+1 

2BAE7230 

CLA 

WORK 

2BAE7240 

ADD 

WORK+1 

2BAE7250 

STO 

WORK 

26AE7260 

CLA 

E-2.2 

2BAE7270 

TSX 

SSQRT.4 

2BAE7280 

ADD 

OCTF 

2BAE7290 

TSX 

FIX. 4 

2BAE7300 

STO 

WORK+1 

2BAE7310 

CAL 

WORK+1 

2BAE7320 

ANA 

MASKR 

2BAE7330 

SLW 

WORK+1 

2BAE7340 

CLA 

WORK 

2BAE7350 

ADD 

WORK+1 

2BAE7360 

STO 

STORGE.l 

STO 

PACKED  WORD 

CONTAINING  3 WORDS 

2BAE7370 

TXI 

*+1.1.1 

2BAE7380 

TXI 

*+1.2.3 

2BAE7390 

TXL 

REP. 2. 119 

2BAE7400 

TRA 

RETURN 

2BAE7410 

RC00E1 

CLA 

MDIR 

REPACK 

518  ROUTINE  FROM  864  FIELD 

2BAE  7420 

SUB 

FID 

2BAE7430 

XCA 

2BAE7440 

MPY 

F10A 

2BAE7450 

XCA 

2BAE7460 

ADD 

MFREO 

2BAE7470 

SUB 

FID 

2BAE7480 

STZ 

WORK 

CHOOSE 

MASK  ROUTINE 

2BAE  7490 

STD 

WORK 

2BAE7500 

CLM 

2BAE7510 

LDO 

WORK 

2BAE7520 

LLS 

0 

2BAE7530 

DVP 

OCT3D 

REMAINDER  IN  AC 

DECR 

2BAE7540 

STO 

WORK 

2BAE7550 
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CLH 

2BAE7560 

LLS 

18 

2BAE7570 

XCA 

2BAE7580 

STD 

RBASE1 

2BAE7S90 

CLA 

WORK 

2BAE7600 

ARS 

18 

2BAE7610 

TZE 

RMASKL 

NO  REMAINDER,  CHOOSE  LEFT  WORD 

2BAE7620 

SUB 

OCTl  A 

2BAE7630 

TN2 

RMASKR 

TRA  IF  REMS,  CHOOSE  RIGHT  WORD 

2BAE  7640 

CLA 

AD12 

REMAINDERS . CHOOSE  CENTER 

2BAE7650 

STA 

RSFTPL 

2BAE  7660 

CLA 

ADMSKC 

2BAE  1670 

STA 

RMSKPL 

2BAE7680 

TRA 

THENPK 

2BAE7690 

RMASKR  CLA 

ADO 

2BAE7700 

STA 

RSFTPL 

2BAE7710 

CLA 

ADMSKR 

2BAE7720 

STA 

RMSKPL 

2BAE7730 

TRA 

THENPK 

2BAE7740 

RMASKL  CLA 

AD24 

2BAE7760 

STA 

RSFTPL 

2BAE  7760 

CLA 

ADMSKL 

2BAE7770 

STA 

RMSKPL 

2BAE7780 

THENPK  A XT 

0,1 

FOR  INPUT  IR1  TO  864 

2BAE7790 

AXT 

0,2 

FOR  OUTPUT  IR2  TO  20720 

2BAfc 7 800 

RBASE1  TXI 

*+1,2,0 

DECR  FROM  RCODE+13 

2BAE7810 

BIGPK  SXD 

INDEX.l 

2BAE7820 

TSX 

SLANOSE ,4 

2BAE7830 

TSX 

INDEX, 0 

2BAE7840 

TNZ 

BGSEPK 

SEA  IF  1 IN  AC 

2BAE7850 

TRA 

RTEST+1 

LAND  POINT 

2BAE7860 

BGSEPK  CLA 

SPACE, 1 

SEA  POINT 

2BAE7870 

TSX 

SS0RT.4 

2BAE7880 

ADD 

OCTF 

2BAE7890 

TSX 

FIX, 4 

2BAE7900 

STO 

WORK 

2BAE7910 

CAL 

WORK 

2BAE7920 

ANA 

MASKR 

2BAE7930 

RSFTPL  ALS 

0 

ADDR  FROM  RMASKR+ 1 ETC 

2BAE7940 

SLW 

WORK 

2BAE7950 

RMSKPL  CAL 

0 

AODR  FROM  RMASKR+3  ETC 

2BAE7960 

COM 

COMPLEMENT  OF  MASK 

2BAE7970 

ana 

STORGE ,2 

clear  place  FOR  PACKED  WORD 

2BAE7980 

ACL 

WORK 

2BAE7990 

SLW 

STORGE, 2 

2BAE8000 

RTFST  TXF 

*+1,2,40 

2BAE8010 

TXI 

*+1.1,1 

2BAE8020 

TXL 

BIGPK, 1.863 

2BAE8030 

TRA 

RETURN 

2BAE8040 

MFREO  HTR 

* 

2BAE8030 

MDIR  HTR 

• 

2BAE8060 

WORK  BSS 

2 

2BAE8070 

OCT3D  OCT 

+000003000000 

2BAE8080 

OCT1A 

+00000  1 

2BAE8090 

t 
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MASKL 

OCT 

-377700000000 

2BAE8100 

MASKC 

OCT 

+000077770000 

2BAE8110 

MASKR 

OCT 

+000000007777 

2BAE8120 

ADMSKL 

HTR 

MASKL 

2BAE8130 

ADMSKC 

HTR 

MASKC 

2BAE8140 

ADMSKR 

HTR 

MASKR 

2BAE8150 

INDEX 

HTR 

0 

2BAE8160 

AD  12 

HTR 

12 

2BAE8170 

ADO 

HTR 

0 

2BAE8180 

AD24 

HTR 

24 

2BAE8190 

OCTF 

OCT 

007000000000 

TO  MODIFY  FLOATING  WORD, BINARY  PT  5-62BAE8200 

F40A 

PZE 

AO 

2BAE821G 

HD 

PZF 

0.0.1 

2BAE8220 

F10A 

PZF 

10 

2BAE8230 

FLOAT 

ORA 

Cl 

NBR  TO  BE  FLOATED  MUST  BE  IN  AC 

ADDR  2BAE8240 

FAD 

Cl 

FLOATED  NBR  STAYS  IN  AC 

2 BAE8250 

TRA 

1 .4 

RETURN  FOR  FLOAT 

2BAE8260 

FIX 

UFA 

Cl 

NBR  TO  BE  FIXED  MUST  BE  IN  AC 

2BAE8270 

LRS 

0 

2BAE8280 

ANA 

Cl  + 1 

2BAE8290 

LLS 

0 

FIXED  NBR  IN  AC  ADDR 

2BAE8300 

TRA 

1 .4 

RETURN  FOR  FIX 

2BAE8310 

Cl 

OCT 

233000000000 

2BAE8320 

OCT 

77777 

2BAE8330 

LSH1FT 

CLA* 

1 .4 

MWHOLE  (0  TO  23) 

2BAE8340 

ARS 

18 

2BAE8350 

ADD 

LCNST 

2BAE8360 

STA 

SI 

ADDRESS  OF  WORD  TO  BE  CHOSEN 

2BAE8370 

CLA* 

2.4 

MP ART  (0  TO  35) 

2BAE8380 

ARS 

18 

2BAE8390 

STA 

S2 

POSITION  IN  WORD 

2BAE8400 

SI 

CAL 

0 

ADDR  FROM  *-4 

2BAE84 10 

S2 

ALS 

0 

ADDR  FROM  *-2 

2BAE8420 

PBT 

2BAE8430 

TRA 

S3 

0 .LAND , PUT  ZERO  IN  AC 

2BAE8440 

CLA 

10CTAD 

1. SEA, PUT  1 IN  AC 

2BAE8450 

TOV 

*+l 

2BAE8460 

TRA 

3.4 

EXIT  FOR  SEA 

2BAE8470 

S3 

CLA 

OOCTAD 

2BAE6480 

TOV 

*+l 

2BAE8490 

TRA 

3.4 

EXIT  FOR  LAND 

2BAE8500 

LCNST 

HTR 

LSTABL 

2BAE8510 

1 OCT  AD 

OCT 

000001000000 

2BAE8520 

OOCTAD 

HTR 

0 

2BAE8530 

LSTABL 

OCT 

+00000000000 

BIT  TABLE 

2BAE8540 

OCT 

+000301034000 

2BAE8550 

OCT 

+006071740000 

BINARY  ONES  REPRESENT  SEA 

2BAE8560 

OCT 

♦143770400003 

ZEROS  FOR  LAND 

2BAE8570 

OCT 

-077600000037 

2BAE8580 

OCT 

-374000000777 

POINTS  ARE  NUMBERED  WEST 

TO  EAST  2BAE8590 

OCT 

-300000017776 

THEN  NORTH  TO  SOUTH.  IE 

SECOND  2BAE8600 

OCT 

+000000377774 

ROW  BEGISN  WITH  NO  33 

2BAE8610 

OCT 

+000007777740 

2BAE8620 

OCT 

+000777774001 

2BAE8630 

B-l8 
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OCT 

-377777700177 

NOTE  SEA  * 5 1# 8 

2BAE8640 

OCT 

-377776037777 

LAND  » 34« £ 

2BAE8650 

OCT 

-377760777777 

TOTAL  ' 864 

2BAE8660 

OCT 

-377417777777 

2BAE8670 

OCT 

-374777777777 

2BAE8680 

OCT 

-317777777776 

2BAE8690 

OCT 

+377777777703 

2BAE8700 

OCT 

-377777774017 

2BAE8710 

OCT 

-377777700377 

2BAE8720 

OCT 

-377776007777 

2BAE8730 

OCT 

-377740177777 

2BAE8740 

OCT 

-377003777777 

2BAE8750 

OCT 

-360077777777 

2BAE8760 

OCT 

-200000000000 

2BAE8770 

SPACE 

COMMON 

864 

2BAE8780 

bigfld 

COMMON 

864 

2BAE8790 

E 

COMMON 

120 

2BAE8800 

STORGF 

COMMON 

20760 

2BAE8810 

EMAX 

COMMON 

120 

2BAE8820 

DIR 

COMMON 

12 

2BAE8830 

FRFO 

COMMON 

10 

2BAE8840 

dfltfo 

COMMON 

10 

2BAE8850 

TABCOL 

COMMON 

25 

2BAE8860 

TABROW 

COMMON 

11 

2BAE8870 

TABLE 

COMMON 

275 

2BAE8880 

JUMP 

COMMON 

864 

2BAE8890 

TIME 

COMMON 

864 

28AE8900 

MI  NO  IS 

COMMON 

12 

2BAE8910 

NEVNTB 

COMMON 

12 

2BAE8920 

NODTAB 

COMMON 

12 

2BAE8930 

END 

2BAE8940 
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Table  C-l  Data:  Case  I (Continued) 


fUEk  CASE  I 
TIME  I 
Rl GMT  SECTION 


-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-c 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

• »• 

-0 

-0 

-0 

-0 

-G 

-0 

-0 

-0 

-0 

0 

-0 

-0 

-C 

-0 

0 

0 

0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

0 

-0 

-0 

-0 

-0 

0 

G 

0 

-0 

-0 

-0 

-C 

-0 

-0 

-0 

0 

0 

0 

-0 

-0 

0 

0 

0 

0 

0 

-0 

-0 

-0 

-0 

-0 

-0 

0 

0 

0 

-0 

-0 

0 

0 

0 

0 

0 

-0 

-0 

-0 

-0 

-0 

0 

0 

0 

0 

0 

0 

0 

0 

-0 

-c 

-0 

0 

-0 

-0 

-0 

-0 

0 

0 

0 

0 

0 

0 

0 

0 

-0 

-c 

-0 

0 

-0 

1? 

-0 

-0 

0 

0 

0 

0 

0 

0 

0 

0 

-0 

-0 

-0 

-0 

-0 

240 

-0 

-0 

0 

0 

0 

0 

0 

0 

0 

0 

-0 

-0 

-0 

-0 

-0 

21 

1 7 

IS 

13 

1 1 

13 

12 

10 

0 

0 

0 

-0 

-0 

-0 

-0 

-0 

2*i  0 

210 

210 

240 

300 

300 

300 

330 

0 

0 

0 

-0 

-0 

-0 

-0 

-0 

25 

20 

17 

17 

17 

17 

17 

16 

14 

12 

0 

-0 

-0 

-0 

-0 

-0 

240 

240 

240 

240 

270 

300 

330 

330 

330 

360 

0 

-0 

-0 

-0 

-0 

-0 

1 rt 

i r 

17 

Id 

19 

20 

21 

22 

21 

2d 

14 

-0 

-0 

-0 

-0 

-0 

240 

240 

240 

240 

270 

270 

300 

3 30 

330 

3 30 

30 

-0 

-0 

-0 

-0 

-0 

13 

U 

IS 

Id 

20 

21 

24 

26 

30 

38 

38 

0 

IS 

IS 

-0 

-0 

24  0 

240 

240 

270 

2 70 

270 

330 

330 

330 

330 

330 

0 

ISO 

ISO 

-0 

-0 

0 

1 1 

14 

Id 

20 

20 

72 

24 

26 

38 

42 

IS 

IS 

16 

16 

-0 

0 

240 

240 

270 

270 

270 

300 

3 30 

330 

330 

330 

2 70 

180 

IdO 

180 

-0 

0 

10 

13 

IS 

1 7 

Id 

19 

20 

24 

30 

42 

32 

Id 

-0 

-0 

-0 

0 

210 

240 

270 

270 

300 

300 

330 

330 

330 

330 

300 

210 

-0 

-0 

-0 

1C 

1 1 

10 

10 

13 

16 

17 

16 

20 

24 

31 

33 

19 

-0 

-0 

-0 

210 

210 

240 

270 

300 

300 

330 

330 

330 

330 

300 

300 

240 

-0 

-0 

'0 

10 

0 

0 

0 

0 

12 

15 

17 

19 

20 

20 

20 

18 

-0 

-0 

'0 

IdO 

0 

0 

0 

0 

360 

360 

360 

330 

330 

300 

270 

240 

-0 

-0 

-0 

0 

0 

0 

0 

0 

0 

13 

16 

17 

18 

18 

18 

17 

IS 

-0 

-0 

0 

0 

0 

0 

0 

0 

360 

360 

360 

330 

300 

270 

270 

240 

-0 

-0 

0 

G 

0 

0 

0 

0 

10 

14 

16 

17 

17 

16 

IS 

14 

-0 

-0 

0 

0 

0 

0 

0 

0 

360 

360 

360 

330 

330 

300 

270 

240 

-0 

-0 

0 

c 

0 

0 

0 

0 

3 

1? 

IS 

16 

16 

IS 

14 

13 

12 

-0 

c 

0 

0 

0 

0 

0 

0 

360 

360 

360 

330 

300 

300 

270 

240 

-0 

0 

0 

3 

0 

0 

3 

0 

*0 

13 

IS 

IS 

14 

IS 

12 

12 

-0 

3 

0 

0 

0 

0 

0 

0 

360 

360 

360 

360 

360 

330 

300 

270 

-0 

0 

0 

0 

0 

0 

0 

0 

0 

1 1 

14 

13 

13 

13 

13 

12 

*0 

0 

0 

0 

0 

c 

0 

0 

0 

30 

360 

360 

360 

36*1 

330 

300 

-0 

0 

0 

0 

0 

0 

0 

0 

0 

11 

13 

14 

14 

14 

14 

-0 

-0 

G 

0 

0 

0 

0 

0 

0 

0 

30 

30 

30 

30 

360 

360 

-0 

-0 

0 

0 

0 

0 

0 

0 

0 

10 

13 

IS 

IS 

IS 

IS 

-0 

-0 

-0 

0 

0 

0 

0 

3 

0 

0 

30 

30 

30 

30 

30 

30 

-0 

-0 

-0 

0 

0 

0 

0 

0 

0 

10 

13 

IS 

IS 

IS 

IS 

IS 

-0 

-0 

-0 

3 

0 

0 

0 

0 

0 

60 

30 

30 

30 

30 

30 

30 

-0 

-0 

-0 

0 

0 

0 

10 

1? 

13 

IS 

16 

16 

16 

16 

IS 

IS 

-0 

-0 

-0 

0 

0 

0 

60 

63 

60 

60 

60 

30 

30 

30 

30 

30 

-0 

-0 

-0 

10 

12 

14 

IS 

16 

16 

16 

16 

16 

16 

IS 

IS 

IS 

-0 

-0 

-0 

90 

90 

60 

60 

60 

60 

60 

60 

60 

30 

30 

30 

30 

-0 

-0 

-0 

13 

IS 

Is 

16 

16 

16 

16 

16 

16 

16 

IS 

14 

12 

-0 

-0 

-0 

90 

9C 

60 

60 

60 

60 

60 

60 

60 

60 

30 

.30 

30 

-0 

-0 

-0 

IS 

16 

16 

16 

17 

17 

17 

17 

16 

IS 

14 

12 

10 

-0 

-0 

-0 

90 

9C 

60 

60 

60 

60 

60 

60 

60 

60 

60 

30 

30 

-0 

-0 

-0 

16 

16 

16 

17 

1 7 

17 

16 

16 

IS 

14 

12 

10 

0 

0 

-0 

-0 

90 

60 

oOO 

60 

60 

60 

60 

60 

60 

60 

60 

600 

0 

0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 
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Table  C-l  Data:  Case  I (Continued) 
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Table  C-l  Data:  Case  I (Continued) 


BAER  CASE  I 
TIME  12 
RIGHT  SECT  I Oft 
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Table  C-l  Data:  Case  I (Continued) 
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Table  C-l  Data;  Case  I (Continued) 


HAER  CASE  I 
TIME  21 
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Table  C-l  Data:  Case  I (Continued) 


•«EI  CASE 
UNC  2* 

LEFT  SECTION 
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Table  C-l  Data:  Case  I (Continued) 

BACK  CASE 
TIME  2% 

RIGHT  SECTION 
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Table  D-l  Data:  Case  II  (Continued) 

8AER  CASE  2 
TIKE  0 
RIGHT  SECT IOM 


-0 

-0 

-0 

-0 

-c 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 
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-0 

-0 

-0 

-0 

0 

0 

0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-c 

65 

-0 

-0 

-0 

-0 

0 

0 

0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

17 

22 

23 

-0 

-0 

0 

0 

0 

0 

15 

-0 

-0 

-0 

-0 

-0 

-0 

60 

65 

65 

-0 

-0 

0 

0 

0 

0 

180 

-c 

-0 

-0 

-0 

-0 

20 

22 

25 

22 

18 

11 

0 

0 

-0 

-0 

-0 

20 

-0 

-0 

-0 

-0 

60 

315 

10 

350 

360 

270 

0 

0 

-0 

-0 

-0 

150 

-0 

SI 

-0 

-0 

29 

26 

23 

21 

16 

11 

0 
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0 

0 

-0 

-0 

-0 

*0 
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28 
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70 
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-0 
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-0 
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180 
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-0 

-0 

-0 

-0 
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28 
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27 

22 

19 

18 

17 

17 
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-0 

-0 

-0 

-0 
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320 

330 

360 

20 

250 

190 

235 

260 

270 

270 

280 

-0 

-0 

-0 

-0 

-0 

29 

36 

35 

36 

36 

30 

25 

23 

20 

19 

19 

19 

21 

15 

-0 

-0 

330 

330 

320 

315 

260 

210 

225 

260 

290 

280 

260 

260 

280 

280 

-0 

-0 

31 

37 

62 

63 

66 

65 

65 

66 

29 

23 

23 

23 

25 

10 

1C 

-0 

330 

320 

315 

270 

270 

260 

220 

260 

270 

270 

270 
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Table  D-l  Data:  Case  II  (Continued) 


mci  case  2 

TIME  1 
RIGHT  SECTION 


-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 
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-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 
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-0 

-0 

-0 
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0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

80 

-0 

-0 

-0 

-0 

0 

0 

0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

19 

2* 

23 

-0 

-0 

15 

10 

0 
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0 

-0 

-0 

-0 

-0 

-0 

-0 

60 

80 

80 

-0 

-0 

130 

150 

0 

0 

0 

-c 

-0 

-0 

-0 

-0 

12 

17 

32 

28 

23 

21 

15 

0 

-0 

-0 

-0 

10 

-0 

-0 

-0 

-0 

60 

60 

80 

80 

100 

130 

160 

0 

-0 

-0 

-0 

60 

-0 

0 

-0 

-0 

18 

2V 

32 

28 

2* 

21 

16 

10 

-0 

-0 

-0 

-0 

-0 

0 

-0 

-0 

60 

80 

90 

100 

no 

150 

160 

180 

-0 

-0 

-0 

-0 

-0 

16 

17 

16 

22 

31 

30 

28 

27 

23 

21 

16 

-0 

-0 

-0 

-0 

-0 

330 

360 
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60 

60 

60 

90 
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180 
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-0 
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-0 
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28 

29 
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32 
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23 
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38 
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31 

20 
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330 
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315 

*00 

290 

270 

260 

225 

210 

180 

-0 

-0 

-0 

-0 

-0 

26 

38 

*8 

52 

S3 
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19 

20 

21 

21 

18 

25 

-0 

-0 

310 
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300 

280 

270 
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220 
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250 

180 

225 

270 
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-0 

28 
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52 

52 

51 

52 

50 

*5 

*0 

33 

28 

23 

18 

25 

25 
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300 
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Table  D-l  Data:  Case  II  (Continued) 


•AM  CASE  2 
TIME  6 

MIGHT  SECTION 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

10 

-0 

-0 

-0 

-0 

IS 

0 

0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

30 

-0 

-0 

-0 

-0 

90 
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0 

-0 

-0 

-0 

' 

-0 

-0 

-0 

-0 

27 

22 

15 

-0 

-0 

0 

12 

0 

0 

0 

-0 

•0 

-0 

-0 

-0 

-0 

360 

10 

30 

-0 

-0 

0 

135 

0 

0 

0 

-c 

-C 

-0 

-0 

20 

35 

35 

25 

17 

10 

0 

13 

-0 

-0 

-0 

0 

.0 

* -o 

-0 

-0 

60 

360 

360 

30 

30 

90 

0 

135 

-0 

-0 

-0 

0 

-0 

0 

-0 

-0 

20 

AO 

A7 

38 

25 

17 

17 

20 

-0 

-0 

-0 

-0 

-0 

■ 

-0 

-0 

360 

360 

360 

360 

360 

220 

160 

200 

-0 

-0 

-0 

-0 

-0 

0 

0 

0 

12 

35 

AS 

52 

AO 

30 

28 

30 

-0 

-0 

-0 

-0 

-0 
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360 
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A2 
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320 
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0 
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0 

0 

-0 

_Q 

c 

0 

0 

0 

0 

0 

0 

0 

0 

70 

60 

60 

0 

0 

-0 

-0 

-c 

-0 

-0 

-0 
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-0 
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Table  D-l  Data:  Case  II  (Continued) 


MW  CASE  2 
TIMC  9 
8I8MT  SECTION 


-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

—0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

27 

-0 

-0 

-0 

-0 

IS 

15 

20 

-0 

-0 
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Table  D-l  Data:  Case  II  (Continued) 

MM  CAM  2 
VIM  IS 
LHV  SCC  T lOM 


-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-c 

-0 

-0 

-0 

-• 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 
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-0 
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-0 
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-0 
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-0 
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-0 

-0 
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-0 

28 
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-0 

-0 

-0 
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Table  D-l  Data:  Case  II  (Continued) 


UH  CASE  2 
TIME  IS 
RI6NT  SECT  1011 


-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-a 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-c 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

0 

-0 

-0 

-0 

-0 

so 

50 

15 

-0 

-c 

-0 

—0 

-0 

-0 

-0 

-0 

0 

-0 

-0 

-0 

-0 

70 

65 

160 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

0 

0 

10 

-0 

-0 

20 

20 

0 

15 

10 

-c 

-0 

-0 

-0 

-0 

-0 

0 

0 

560 

-0 

-0 

550 

65 

0 

160 

155 

-0 

-0 

-0 

-0 

-0 

0 

0 

0 

27 

16 

25 

65 

62 

-0 

-0 

-0 

25 

-0 

-0 

-0 

-0 

0 

0 

0 

560 

560 

560 

550 

560 

-0 

-0 

-0 

190 

-0 

25 

-0 

-0 

20 

20 

21 

60 

50 

62 

55 

60 

-0 

-0 

-0 

-c 

-0 

1 95 

-0 

-0 

210 

270 

500 

550 

560 

560 

500 

520 

-0 

-0 

-0 

-0 

-0 

AO 

62 

65 

65 

57 

15 

55 

65 

60 

55 

60 

-0 

-0 

-0 

-0 

-0 

60 

160 

160 

200 

270 

500 

520 

550 

560 

510 

500 

-0 

-0 

-0 

-0 

-0 

55 

55 

50 

79 

2V 

29 

25 

66 

65 

52 

58 

-0 

-0 

-0 

-0 

-0 

155 

170 

150 

150 

225 

500 

515 

520 

520 

510 

270 

-0 

-0 

-0 

-0 

-0 

2* 

27 

27 

25 

21 

20 

15 

29 

59 

65 

55 

-0 

-0 

-0 

-0 

-0 

2 TO 

200 

160 

160 

100 

270 

515 

520 

515 

500 

2T0 

-0 

-0 

-0 

-0 

-0 

10 

25 

51 

51 

20 

15 

12 

21 

50 

62 

62 

52 

25 

50 

-0 

-0 

270 

25C 

160 

160 

160 

150 

515 

515 

515 

500 

290 

250 

220 

260 

-c 

-0 

20 

IS 

65 

55 

25 

10 

0 

21 

55 

59 

62 

55 

25 

60 

6C 

-0 

270 

260 

160 

160 

160 

150 

0 

515 

515 

500 

290 

2 70 

220 

250 

250 

-0 

25 

10 

55 

55 

27 

20 

10 

12 

20 

52 

50 

51 

50 

-0 

-0 

-0 

260 

250 

150 

160 

150 

150 

100 

515 

515 

500 

290 

290 

250 

-0 

-0 

-0 

20 

22 

51 

55 

50 

21 

15 

25 

51 

20 

20 

25 

50 

-0 

-0 

-0 

270 

225 

225 

160 

150 

150 

160 

150 

515 

500 

270 

290 

250 

-0 

-0 

-0 

25 

10 

50 

55 

27 

19 

11 

IS 

19 

15 

15 

18 

20 

-0 

-0 

-0 

225 

225 

250 

260 

150 
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160 
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-0 

-0 

•0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 
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-0 

-0 

-0 

-0 

-0 
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-0 

-0 

-0 

d-14 
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Table  D-l  Data:  Case  II  (Continued) 


BAHI  CASE  2 

TIME  18 

IE8T  SECTION 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

«0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

0 

0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

0 

0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

0 

0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

0 

0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

0 

22 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

0 

60 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

25 

31 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

10 

90 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

3S 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

70 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

35 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

180 

-0 

-0 

-0 

-0 
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Table  D-l  Data:  Case  II  (Continued) 
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APPENDIX  E 


NEW  YORK  UNIVERSITY 

COLLEGE  OF  ENGINEERING 

UNIVERSITY  HEIGHTS.  NEW  YORK  53.  N Y. 


DEPARTMENT  OF  METEOROLOGY 
AND  OCEANOGRAPHY 


3 November  1961 


TELEPHONE:  LUDLOW  4-0700 


'■*r.  Ledolph  Baer, 

10?0  Havre  Court 
Sunnyvale,  Calif, 

Dear  Mr.  Baert 

The  purpose  of  this  letter  le  to  describe  to  you  the  nay  In 
which  the  data  that  we  have  furnished  you  with  reference  to  the  storm 
ever  the  *’orth  Atlantic  in  December  1959  were  obtained.  The  collection 
of  this  data  has  been  a cooperative  effort  of  a number  of  people  who 
are  Interested  in  the  study  of  gravity  waves.  They  are 

Dr.  C.L.  Bretachnelder 

Dr.  J.  Darbyshlre 

Dr.  H.L.  Crutcher 

Dr.  H.  Klche 

Dr.  H.  Walden 

Dr.  B.V.  Wilson 

Dr.  W.J. Pierson,  Jr. 

Dr.  G.  Neumann 

The  original  data  were  supplied  to  me  by  the  National  In- 
stitute of  Oceanography.  They  were  digitised  with  funds  made  available 
by  the  Office  at  Naval  Research  on  their  contract  Nonr  ?85(03).  The 
qpectra  were  ccegroted  at  New  Tork  University  by  Bnanuel  Mehr.  Supplement- 
ary wind  data  were  searched  at  the  National  Weather  Records  Center  by 
Dr.  Harold  Crutcher  and  provided  to  us  to  supplement  the  wind  fields 
that  Were  deduced  from  both  the  British  Weather  Maps  and  from  our 
Northern  Hemisphere  6 hourly  surface  charts. 

It  Is  the  intent  of  the  above  named  group  to  publish  these 
data  as  a report  on  a particularly  severe  storm  In  the  North  Atlantic 
ocean.  Just  when  this  publication  will  appear  is  not  yet  known  ae  the 
paper  has  not  been  written,  but  all  of  the  data  for  the  paper  Is  now 
available.  Perhaps  at  this  time  each  of  the  Interested  parties  will 
then  prepare  forecasts  of  the  waves  based  on  the  techniques  with  which 
he  Is  most  familiar  to  see  how  well  the  various  methods  work  out. 


WJPjEM 


Sincerely  yours. 


Willard  J.  Pierson,  Jr. 
Professor  of  Meteorology 
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